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1. THE "MINIMUM" PRINCIPLES 



Forward: The history of "minimum" principles in physics is long and inter- 
esting. The study of such principles is based on the idea that the nature acts 
always in such a way that the important physical quantities are minimized 
whenever a real physical process takes place. The mathematical background 
for these principles is the variational calculus. 
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1. Introduction 

The empirical evidence has shown that the motion of a particle in an inertial 
system is correctly described by Newton's second law F = dp/dt, whenever 
possible to neglect the relativistic effects. When the particle happens not to 
be forced to a complicated motion, the Cartesian coordinates are sufficient 
to describe the movement. If none of these conditions are fulfilled, rather 
complicated equations of motion are to be expected. 

In addition, when the particle moves on a given surface, certain forces called 
constraint forces must exist to maintain the particle in contact with the 
surface. Such forces are not so obvious from the phenomenological point 
of view; they require a separate postulate in Newtonian mechanics, the one 
of action and reaction. Moreover, other formalisms that may look more 
general have been developed. These formalisms are equivalent to Newton's 
laws when applied to simple practical problems, but they provide a general 
approach for more complicated problems. The Hamilton principle is one 
of these methods and its corresponding equations of motion are called the 
Euler-Lagrange equations. 

If the Euler-Lagrange equations are to be a consistent and correct descrip- 
tion of the dynamics of particles, they should be equivalent to Newton's 
equations. However, Hamilton's principle can also be applied to phenom- 
ena generally not related to Newton's equations. Thus, although HP does 
not give a new theory, it unifies many different theories which appear as 
consequences of a simple fundamental postulate. 



The first "minimum" principle was developed in the field of optics by Heron 
of Alexandria about 2,000 years ago. He determined that the law of the 
reflection of light on a plane mirror was such that the path taken by a light 
ray to go from a given initial point to a given final point is always the shortest 
one. However, Heron's minimum path principle does not give the right law 
of reflection. In 1657, Fermat gave another formulation of the principle by 
stating that the light ray travels on paths that require the shortest time. 
Format's principle of minimal time led to the right laws of reflection and 
refraction. The investigations of the minimum principles went on, and in 
the last half of the XVII century, Newton, Leibniz and Bernoulli brothers 
initiated the development of the variational calculus. In the following years, 
Lagrange (1760) was able to give a solid mathematical base to this principle. 
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In 1828, Gauss developed a method of studying Mechanics by means of his 
principle of minimum constraint. Finally, in a sequence of works published 
during 1834-1835, Hamilton presented the dynamical principle of minimum 
action. This principle has always been the base of all Mechanics and also of 
a big part of Physics. 

Action is a quantity of dimensions of length multiplied by the momentum 
or energy multiplied by time. 

2. The action principle 

The most general formulation of the law of motion of mechanical systems 
is the action or Hamilton principle. According to this principle every me- 
chanical system is characterized by a function defined as: 

or shortly L ^q, q, tj , and the motion of the system satisfies the following 

condition: assume that at the moments ti and t2 the system is in the posi- 
tions given by the set of coordinates (j'^^' y (j^^^; the system moves between 
these positions in such a way that the integral 

S = J\(^q,q,t^dt (1) 

takes the minimum possible value. The function L is called the Lagrangian 
of the system, and the integral (1) is known as the action of the system. 

The Lagrange function contains only q and q, and no other higher-order 
derivatives. This is because the mechanical state is completely defined by 

its coordinates and velocities. 

Let us establish now the difererential equations that determine the minimum 
of the integral (1). For simplicity we begin by assuming that the system has 
only one degree of freedom, therefore wc arc looking for only one function 
q (t). Let q = q{t) be the function for which is a minimum. This means 
that S grows when one q (t) is replaced by an arbitrary function 

q{t)+Sq{t), (2) 

where Sq (t) is a small function through the interval from ti to t2 [it is called 
the variation of the function q{t)]. Since at ti and t2 all the functions (2) 
should take the same values q^^^ and q^'^\ one gets: 
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5q (ti) = 5q (ts) = 0. 
What makes S change when q is replaced by q + 5q is given by: 



(3) 



j\{q + 5q,q+5q,t^dt- [q, q, t 



dt. 



An expansion in series of this difference in powers of Sq and S q begins by 
terms of first order. The necessary condition of minimum (or, in general, 
extremum) for S is that the sum of all terms turns to zero; Thus, the action 
principle can be written down as follows: 



SS = S 
or by doing the variation: 

r*i / dL 



t2 



L[q,q,t]dt = 0, 



(4) 



/ 

Jt2 



^ 5q + 6q]dt = 

^1 dq ) 



Taking into account that 5 q= d/dt{dq), we make an integration by parts 
to get: 



SS 



dL 
Idq 



Sq 



t2 



+ 



*i dL 



t2 



d dL\ ^ , ^ 
dq dtgqj 



(5) 



Considering the conditions (3), the first term of this expresion disappears. 
Only the integral remains that should be zero for all values of Sq. This is 
possible only if the integrand is zero, which leads to the equation: 

dL _ d dL _^ 
dq dt Q q 

For more degrees of freedom, the s different functions qi{t) should vary 
independently. Thus, it is obvious that one gets s equations of the form: 



d dL 



dL 
dqi 



= 



(z = l,2,...,s) 



(6) 



dt \dq,, 

These are the equations we were looking for; in Mechanics they axe called 
Euler-Lagrange equations. If the Lagrangian of a given mechanical system 
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is known, then the equations (6) form the relationship between the accel- 
erations, the velocities and the coordinates; in other words, they are the 
equations of motion of the system. From the mathematical point of view, 
the equations (6) form a system of s differential equations of second order 
for s unknown functions qi{t). The general solution of the system contains 
2s arbitrary constants. To determine them means to completely define the 
movement of the mechanical system. In order to achieve this, it is neces- 
sary to know the initial conditions that characterize the state of the system 
at a given moment (for example, the initial values of the coordinates and 
velocities. 

3. D'Alembert principle 
The virtual displacement of a system is the change in its configurational 
space under an arbitrary infinitesimal variation of the coordinates 6ri, which 
is compatible with the forces and constraints imposed on the system at the 
given instant t. It is called virtual in order to distinguish it from the real 
one, which takes place in a time interval dt, during which the forces and the 
constraints can vary. 

The constraints introduce two types of difficulties in solving mechanics prob- 
lems: 

(1) Not all the coordinates are independent. 

(2) In general, the constraint forces are not known a priori; they are some 
unknowns of the problem and they should be obtained from the solution 
looked for. 

In the case of holonomic constraints the difficulty (1) is avoided by introduc- 
ing a set of independent coordinates {qi,q2,...,Qm, where m is the number of 
degrees of freedom involved). This means that if there are m constraint equa- 
tions and SN coordenates (xi, ...,X3n), we can eliminate these n equations 
by introducing the independent variables {qi,q2, ,Qn)- A transformation 
of the following form is used 

xi = fi{qi,:;qm,t) 



X3N = f3N{qi, qn, t) , 

where n = 3N — m. 

To avoid the difficulty (2) Mechanics needs to be formulated in such a way 
that the forces of constraint do not occur in the solution of the problem. 
This is the essence of the '^principle of virtual worM^ . 
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Virtual work: We assume that a system of N particles is described by 
3N coordenates {xi,X2, ...,X3n) and let -Fi,i^2,...,-^37V be the components of 
the forces acting on each particle. If the particles of the system display 
infinitesimal and instantaneous displacements 6x1,6x2, ■■■ySx^isi under the 
action of the 3N forces, then the performed work is: 

3N 

SW = J2Fj6xj. (7) 

Such displacements are known as virtual displacements and SW is called 
virtual work, {7) can be also written as: 

TV 

6W = Y,Fa-Sr. (8) 

a=l 

Forces of constraint: besides the applied forces Fa \ the particles can be 
acted on by forces of constraint Fq. 

The principle of virtual work: Let Fq, be the force acting on the particle 

(e) 

a of the system. If we separate Fq, in a contribution from the outside Fa 
and the constraint Ra 

Fa = F(f ) + Ra ■ (9) 
and if the system is in equilibrium, then 

F, = F(f) + = . (10) 
Thus, the virtual work due to all possible forces Fq, is: 

N N 

W=J2Fa-Sra = Yl (^i'^ + Ra) ■ SVa = . (11) 
a=l a=l 

If the system is such that the constraint forces do not make virtual work, 
then from (11) we obtain: 

N 

Y: Fi«) ■Sra = 0. (12) 

a=l 

Taking into account the previous definition, we are now ready to introduce 
the D'Alembert principle. According to Newton, the equation of motion is: 
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and can be written in the form 



Fa- Pa= , 



which tehs that the particles of the system would be in equilibrium under the 

action of a force equal to the real one plus an inverted force — Pj. Instead 
of (12) we can write 



a=i ^ ' 

and by doing the same decomposition in applied and constraint forces (fa), 
we obtain: 



Again, let us limit ourselves to systems for which the virtual work due to 
the forces of constraint is zero leading to 



which is the D'Alembert's principle. However, this equation does not have 
a useful form yet for getting the equations of motion of the system. There- 
fore, we should change the principle to an expression entailing the virtual 
displacements of the generalized coordinates, which being independent from 
each other, imply zero coefficients for Sq^a- Thus, the velocity in terms of 
the generalized coordinates reads: 




(13) 




a=l 




(14) 




where 



r, 




t) 



Similarly, the arbitrary virtual displacement 6ra can be related to the virtual 
displacements Sqj through 



Then, the virtual work expressed in terms of the generalized coordinates 
will be: 



N 



^ F„ • 5r„ = ^ F„ • -^Sqj = ^ QjSqj , 



(15) 



a=l 



where the Qj are the so-called components of the generalized force, defined 
in the form 

Now if we see eq. (14) as: 

P -dVa = X! "^a Ta -dVa (16) 
a a 

and by substituting in the previous results we can see that (16) can be 
written: 



d 



dVa 



E 



d dT 



dt \dq 



Qj. 



3-' 

(17) 

The variables qj can be an arbitrary system of coordinates describing the 
motion of the system. However, if the constraints are holonomic, it is possi- 
ble to find systems of independent coordinates qj containing implicitly the 
constraint conditions already in the equations of transformation Xi = fi if 
one nullifies the coefficients by separate: 



6qj = , 



dt \QqJ Oqa 



(18) 



There are m equations. The equations (18) are sometimes called the La- 
grange equations, although this terminology is usually applied to the form 
they take when the forces are conservative (derived from a scalar potential 
V) 



Fa = -v^y. 



Then Qj can be written as: 
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The equations (18) can also be written in the form: 



d_ (dT_\ _ d{T - V) 



(19) 



and defining the Lagrangian L in the form L = T — V one gets 




d_ / dL 



) 




(20) 



These are the Lagrange equations of motion. 

4. - Phase space 

In the geometrical interpretation of mechanical phenomena the concept of 
phase space is very much in use. It is a space of 2s dimensions whose axes 
of coordinates are the s generalized coordinates and the s momenta of the 
given system. Each point of this space corresponds to a definite mechanical 
state of the system. When the system is in motion, the representative point 
in phase space performs a curve called phase trajectory. 

5. - Space of configurations 

The state of a system composed of n particles under the action of m con- 
straints connecting some of the 3n cartesian coordinates is completely de- 
termined by s = 3n — m generalized coordinates. Thus, it is possible to 
descibe the state of such a system by a point in the s dimensional space 
usually called the configuration space, for which each of its dimensions cor- 
responds to one qj. The time evolution of the system will be represented by a 
curve in the configuration space made of points describing the instantaneous 
configuration of the system. 

6. - Constraints 

One should take into account the constraints that act on the motion of the 
system. The constraints can be classified in various ways. In the general 
case in which the constraint equations can be written in the form: 
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where the Cai are functions of only the coordinates (the index a counts 
the constraint equations). If the first members of these equations are not 
total derivatives with respect to the time they cannot be integrated. In other 
words, they cannot be reduced to relationships between only the coordinates, 
that might be used to express the position by less coordinates, corresponding 
to the real number of degrees of freedom. Such constraints axe called non 
holonomic (in contrast to the previous ones which are holonomic and which 
connect only the coordinates of the system ). 

7. Hamilton's equations of motion 
The formulation of the laws of Mechanics by means of the Lagrangian as- 
sumes that the mechanical state of the system is determined by its gener- 
alized coordinates and velocities. However, this is not the unique possible 
method; the equivalent description in terms of its generalized coordinates 
and momenta has a number of advantages. 

Turning from one set of independent variables to another one can be achieved 

by what in mathematics is called Legendre transformation. In this case the 
transformation takes the following form where the total differential of the 
Lagrangian as a function of coordinates and velocities is: 

■sr^ dL , y-^ dL ■ 

i i dqi 

that can be written as: 



dL = Ypidqi + Ypidqi , (21) 

i i 

where we already know that the derivatives dL/d qi, are by definition the 

generalized momenta and moreover dL / dqi =pi by Lagrange equations. 
The second term in eq. (21)can be written as follows 



By attaching the total diferential d {^J^Pi lij to the first term and changing 
the signs one gets from (21): 

d (j2Pi '^i = - E Pi dqi + Y,Pi 1i ■ (22) 
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The quantity under the diferential is the energy of the system as a func- 
tion of the coordinates and momenta and is called Hamiltonian function or 
Hamiltonian of the system: 

H{jp,q,t) = Y,Pi1i-L . (23) 

i 

Then from ec. (22) 

dH = -Y,P, dqi+J2pi Qi 

where the independent variables are the coordinates and the momenta, one 
gets the equations 

9H ■ dH 

dpi dqi 
These are the equations of motion in the variables q y p and they are called 
Hamilton's equations. 

8. Conservation laws 
8.1 Energy 

Consider first the conservation theorem resulting from the homogeneity of 
time. Because of this homogeneity, the Lagrangian of a closed system does 
not depend explicitly on time. Then, the total time diferential of the La- 
grangian (not depending explicitly on time) can be written: 

dt ^dqi ^ Q q. 

and according to the Lagrange equations we can rewrite the previous equa- 
tion as follows: 

dL ^ ■ d I dL\ dL ■■ ^ d / ■ dL\ 

E'^^:^ — +E — 



dt ■ dt \dq ) i dq i ^t \ q q 



or 



Prom this one concludes that the quantity 



E^J^Ii — 'L (25) 
i dqi 
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remains constant during the movement of the closed system, that is it is 
an integral of motion. This constant quantity is called the energy E of the 
system. 

8.2 Momentum 

The homogeneity of space implies another conservation theorem. Because 
of this homogeneity, the mechanical properties of a cosed system do not 
vary under a parallel displacement of the system as a whole through space. 
We consider an infinitesimal displacement e (i.e., the position vectors r„ 
turned into + e) and look for the condition for which the Lagrangian does 
not change. The variation of the function L resulting from the infinitesi- 
mal change of the coordinates (maintaining constant the velocities of the 
particles) is given by: 



extending the sum over all the particles of the system. Since e is arbitrary, 
the condition = is equivalent to 



and taking into account the already mentioned equations of Lagrange 
^ d f dL\ d ^ dL „ 



Thus, we reach the conclusion that for a closed mechanical system the vec- 
torial quantity called impetus /momentum 



remains constant during the motion. 
8.3 Angular momentum 

Let us study now the conservation theorem coming out from the isotropy of 
space. For this we consider an infinitesimal rotation of the system and look 
for the condition under which the Lagrangian does not change. 
We shall call an infinitesimal rotation vector 54) a vector of modulus equal 
to the angle of rotation S(f) and whoose direction coincide with that of the 
rotation axis. We shall look first to the increment of the position vector 





(26) 
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of a particle in the system, by taking the origin of coordinates on the axis 
of rotation. The hneal displacement of the position vector as a function of 
angle is 

\5r\ = r sin 95(p , 

(see the figure) . The direction of the vector Sr is perpendicular to the plane 
defined by r and 6(p, and therefore, 

Sr =5^ X r . (27) 



A 



o 




The rotation of the system changes not only the directions of the position 
vectors but also the velocities of the particles that are modified by the same 
rule for all the vectors. The velocity increment with respect to a fixed frame 
system will be: 

6-v = Scf) X Y . 

We apply now to these expressions the condition that the Lagrangian does 
not vary under rotation: 

«^i:(t-*^»+|;-*v.)^o 

and substituting the definitions of the derivatives dLjdWa por Pa and dLjdva 
from the Lagrange equations by we get 

^ (Pa -50 X r„ + p„ • (5^ X = , 

a ^ ' 
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or by circular permutation of the factors and getting 5^ out of the sum: 
because Scf) is arbitrary, one gets 

a 

Thus, the conclusion is that during the motion of a closed system the vec- 
torial quantity called the angular (or kinetic) momentum is conserved. 

M = J2^aXPa ■ 

a 



9.- Applications of the action principle 

a) Equations of motion 

Find the eqs of motion for a pendular mass sustained by a resort, by directly 
applying Hamilton's principle. 




For the pendulum in the figure the Lagrangian function is 
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therefore 
rt2 



1 .2 -2 I 

L = -m(r ) + mgr cos 9 — -k{'i — ro)^ , 



ft2 rt2 

/ 5Ldt = 
Jti Jti 



.2 



m [r S r +r +r 9 S 0] + mgSr cos — mgrS9 sin 9 — A;(r — ro)Sr 



dt 



mrdrdt = mr d{6r) = d (m r drj — m6r r dt 
In the same way 



mr 



, . d \ mr^ 

9dt = d [mr^ 6 6 9j - ^G^-^^ — '-dt 



d mr^ 9 S 9] -S9 ( mr^ 9 +2mr r0 dt 



Therefore, the previous integral can be written 



rt2 
Jti 



.2 



m r —mr —mg cos 9 + k {r — rg) ? + { mr 9 +2mr r0 +mgr sin 



dt 



t2 



d(mr dr^ +d imr'^9'^ 9 50 



= . 



Assuming that both 5r and 59 are equal zero at ti and t2, the second integral 
is obviously nought. Since 5r and 50 are completely independent of each 
other, the first integral can be zero only if 



m r —mr 9 —mg cos + k{r — rg) = 



and 



mr +2mr rQ +mgrsm9 = , 
These are the equations of motion of the system. 
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b) Exemple of calculating a minimum value 

Prove that the shortest Une between two given points pi and p2 on a cihnder 
is a helix. 

The length S of an arbitrary hne on the cihnder between pi and p2 is given 

by 



/ 



P2 
Pi 



az 



1/2 



dz , 



where r, 6 and z are the usual cilindrical coordinates for r = const. A 
relationship between 9 and z can be determined for which the last integral 
has an extremal value by means of 



where </>=[! + r^O'^] ^^'^ y 6' = ^, but since d(t)/de = we have 



Jz yae' 



de 







^^ = (l + r'e''y"\^e' = c,= const., 

therefore rO' = ci- Thus, rQ = ciz + ca, which is the parametric equation of 
a helix. Assuming that in p\ we have ^ = and 2; = 0, then C3 = 0. In pi, 
make 6 = 02 and z = Z2, therefore C2 = r92/z2, and r9 = {r92/z2)z is the 
final equation. 



References 



L. D. Landau and E. M Lifshitz, Mechanics, Theoretical Physics, vol I, 
(Pergammon, 1976) 

H. Goldstein, Classical Mechanics, (Addison- Wesley, 1992) 



18 



2. MOTION IN CENTRAL FORCES 

Forward: Because of astronomical reasons, the motion under the action of 
central forces has been the physical problem on which pioneer researchers 
focused more, either from the observational standpoint or by trying to disen- 
tangle the governing laws of motion. This movement is a basic example for 
many mathematical formalisms. In its relativistic version, Kepler's problem 
is yet an area of much interest. 

CONTENTS: 

2.1 The two-body problem: reduction to the one-body problem 

2.2 Equations of motion 

2.3 Differential equation of the orbit 

2.4 Kepler's problem 

2.5 Dispertion by a center of forces (with example) 
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2.1 Two-body problem: Reduction to the one-body problem 

Consider a system of two material points of masses mi and m2, in which 
there are forces due only to an interaction potential V. We suppose that V 
is a function of any position vector between mi and m2, r2 — ri, or of their 
relative velocities r2 — ri, or of the higher-order derivatives of r2 — ri. Such 
a system has 6 degrees of freedom and therefore 6 independent generalized 
coordinates. 

We suppose that these are the vector coordinates of the center-of-mass R, 
plus the three components of the relative difference vector r = r2 — ri. The 
Lagrangian of the system can be written in these coordinates as follows: 

L = r(R,r)-y(r,r,i=, ). (1) 

The kinetic energy T is the sum of the kinetic energy of the center-of-mass 
plus the kinetic energy of the motion around it, T': 



being 



T= ^(mi + m2)R^ + r, 



T= -miri + -m2r2. 



Here, ri' and r2' are the position vectors of the two particles with respect to 
the center-of-mass, and they are related to r by means of 

m2 , mi 
ri = ■ r, r2 = ■ r . 2) 

mi -|- m2 mi -|- m2 

Then, Intakes the form 

1 mim2 .2 
1 = r 

2 mi -|- m2 

and the total Lagrangian as given by equation (1) is: 

L=\{m,+ m2)R2 + l^!^r2 _ y(r, r, r, ) , (3) 

2 2 mi -|- m2 

where from the reduced mass is defined as 

mim2 .1 1 1 
/i = o — = 1 . 

mi -|- m2 /U mi m2 

Then, the equation (^) can be written as follows 

L = ^(mi + m2)R^ + ^/uf^ - V{y, r, r, ). 
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From this equation we see that the coordinates R, are cycHc implying that 
the center-of-mass is either fixed or in uniform motion. 
Now, none of the equations of motion for r will contain a term where R or 
R will occur. This term is exactly what wc will have if a center of force 
would have been located in the center of mass with an additional particle 
at a distance r away of mass n (the reduced mass). 

Thus, the motion of two particles around their center of mass, which is due 
to a central force can be always reduced to an equivalent problem of a single 
body. 

2.2 Equations of motion 

Now we limit ourselves to conservative central forces for which the potential 
is a function of only r, V{r), so that the force is directed along r. Since 
in order to solve the problem we only need to tackle a particle of mass 
m moving around the fixed center of force, we can put the origin of the 
reference frame there. As the potential depends only on r, the problem has 
spherical symmetry, that is any arbitrary rotation around a fixed axis has no 
effect on the solution. Therefore, an angular coordinate representing that 
rotation should be cyclic providing another considerable simplification to 
the problem. Due to the spherical symmetry, the total angular momentum 

L = r X p 

is conserved. Thus, it can be inferred that r is perpendicular to the fixed 
axis of L. Now, if L = the motion should be along a line passing through 
the center of force, since for L = r and r are parallel. This happens only 
in the rectilinear motion, and therefore central force motions proceed in one 
plane. 

By taking the z axis as the direction of L, the motion will take place in the 
(x, y) plane. The spherical angular coordinate (p will have the constant value 
tt/2 and we can go on as foollows. The conservation of the angular momen- 
tum provides three independent constants of motion. As a matter of fact, 
two of them, expressing the constant direction of the angular momentum, 
are used to reduce the problem of three degrees of freedom to only two. The 
third coordinate corresponds to the conservation of the modulus of L. 
In polar coordinates the Lagrangian is 

L = ^m{f^ + r^e^) - V{r) . (4) 
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As we have seen, ^? is a cyclic coordinate whose canonicaUy conjugate mo- 
mentum is the angular momentum 

Pe = — r = mr u , 
89 

then, one of the equations of motion will be 

p, = |(mr2^) = 0. (5) 

This leads us to 

= I = cte , (6) 

where I is constant modulus of the angular momentum. Prom equation (5) 
one also gets 

The factor 1/2 is introduced because (r^0)/2 is the areolar velocity (the area 
covered by the position vector per unit of time). 

The conservation of the angular momentum is equivalent to saying that the 
areolar velocity is constant. This is nothing else than a proof of Kepler's 
second law of planetary motion: the position vector of a planet covers equal 
areas in equal time intervals. However, we stress that the constancy of the 
areolar velocity is a property valid for any central force not only for inverse 
square ones. 

The other Lagrange equation for the r coordinates reads 

— (mr)-mre^ + — = 0. (8) 
at or 

Denoting the force by /(r), we can write this equation as follows 

mr — mrO"^ = f(r) . (9) 
Using the equation (6), the last equation can be rewritten as 

mf--^ = /(r). (10) 
Recalling now the conservation of the total energy 

E = T + V = ^m{f'^ + r^9^) + V{r) . (11) 
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we say that is a constant of motion. This can be derived from the equa- 
tions of motion. The equation (10) can be written as follows 



mr 



d_ 

dr 



V{r) + - 



1 f 



2 mr^ 



(12) 



and by multiplying by r both sides, we get 



mrr 



d 1 

— (—mr) 
dV2 ' 



d_ 
di 



Vir) + 



1 f 



2 mr^ 



or 



d_ 

dt 



-mr +V{r) + - ;r 

2 2 mr^ 



Thus 



-mr +V[r) + - ^ 

2 2 mr^ 



(13) 



and since (/^/2mr^) = {mr'^9/2), the equation (13) is reduced to (11). 
Now, let us solve the equations of motion for r and 9. Taking r from equation 
(|l3|) , we have 



m 



-{E-V 



P 



2mr^ 



) 



or 



dt 



dr 



2/_2_ 

m 



{E-V 



P - 

2mr^- 



(14) 



(15) 



Let To be the value of r at t 
equation reads 

t -- 



0. The integral of the two terms of the 
dr 



2/2. 
m 



[E-V 



(16) 



This equation gives t as a function of r and of the constants of integration 
E, I and tq. It can be inverted, at least in a formal way, to give r as a 
function of t and of the constants. Once we have r, there is no problem to 
get 6 starting from equation (6), that can be written as follows 



dO 



Idt 



mr^ 



If ^0 is the initial value of 9, then (17) will be 

dt 



I 



mr^{t) 



(17) 



(18) 
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Thus, we have already get the equations of motion for the variables r and 
6. 



2.3 The differential equation of the orbit 

A change of our standpoint regarding the approach of real central force prob- 
lems prove to be convenient. Till now, solving the problem meant seeking r 
and 9 as functions of time and some constants of integration such as E, I, 
etc. However, quite often, what we are really looking for is the equation of 
the orbit, that is the direct dependence between r and 9, by eliminating the 
time parameter t. In the case of central force problems, this elimination is 
particularly simple because t is to be found in the equations of motion only 
in the form of a variable with respect to which the derivatives are performed. 
Indeed, the equation of motion (^) gives us a definite relationship between 
dt and d9 

Idt = mr^d9. (19) 

The corresponding relationship between its derivatives with respect to t and 
9 is 

d I d 

dt^^d9' ^ ' 

This relationship can be used to convert ( p^ ) in a differential equation for 
the orbit. At the same time, one can solve for the equations of motion and 
go on to get the orbit equation. For the time being, we follow up the first 
route. 

From equation ( pOD we can write the second derivative with respect to t 

d^ d I f d I 



dt"^ d9 mr"^ \d9 mr"^ 
and the Lagrange equation for r, (|^), will be 
I d f I dr\ I 



d9 V fnT"^ d9 J mr^ 
But 

l_dT _ d{l/r) 
^d9 ~ d9~ ■ 
Employing the change of variable u = 1/r, we have 



fir) . (21) 




-/ ( - ) • (22) 



u 
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Since 

d dr d Id 
du dO dr v? dr ' 
equation (^) can be written as follows 

d^u m d f 1 



dO^ Pdu 



+ ^ = -n:j^y - ■ (23) 



Any of the equations (|22D or ( |23D is the differential equation of the orbit if 
we know the force / or the potential V. Vice versa, if we know the orbit 
equation we can get / or V. 

For an arbitrary particular force law, the orbit equation can be obtained 
by integrating the equation (|2^). Since a great deal of work has been done 
when solving (|To|), we are left with the task of eliminating t in the solution 
(|l5|) by means of (lir 



Idr , ,, 

d9 = 1^^^=^^^= , (24) 



mr^ ■ ^ 



E-V{r)-^ 



or 



"^"^ +^0. (25) 



2/2mE _ 2mU _ 1 



By the change of variable u = 1/r, 

du 

^0 



no ^2mB _ 2ml/ _ ^2 

which is the formal solution for the orbit equation. 



(26) 



2.4 Kepler's problem: the case of inverse square force 

The inverse square central force law is the most important of all and therefore 
we shall pay more attention to this case. The force and the potential are: 

f = -^2 y y = -^- (27) 

To integrate the orbit equation we put (23) in (22), 
d'^u mf(l/u) mk 
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Now, we perform the change of variable y = u — ^ , in order that the 
differential equation be written as follows 

possessing the solution 

y = Bcos{e-e) , 

where B and ^'are the corresponding integration constants. The solution in 
terms of r is 

l = !I^[l + ecos(0-0)], (29) 

where 

mk 

We can get the orbit equation from the formal solution (p6|). Although the 
procedure is longer than solving the equation (^8|) , it is nevertheless to do it 
since the integration constant e is directly obtained as a function of E and 

We write equation (]26|) as follows 

(30) 



2 2mE 2mU „,2 
1^ P ^ 



where now one deals with a definite integral. Then 0'of ( pOf ) is an integration 
constant determined through the initial conditions and is not necessarily the 
initial angle 6q at t = 0. The solution for this type of integrals is 



/ 



dx 1 

arccos 



where 



l3 + 2jx 



(31) 



q = 0^ - 4a7. 

In order to apply this type of solutions to the equation ( |30| ) we should make 

2mE 2mk 



and the discriminant q will be 



f2mk\ . ^ 



2Ef 



\ P J \ mk"^ 
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With these substitutions, (pO|) is 



arccos 



mk 



^'1 + 



2EP 
mk-^ 



For u = 1/r, the resulting orbit equation is 
1 mk 



2EP 

1 + \H — cos( 



m. 



A;2 



(32) 



Comparing (|3^) with the equation (p^) we notice that the value of e is: 



fc2 



(33) 



The type of orbit depends on the value of e according to the following table: 



e > 1, E >0 

e = l, ^ = 

e < 1, ^ < 

e = = - 



mk^ 



hyperbola, 

parabola, 

elipse, 

circumference. 



2.5 Dispersion by a center of force 

From a historical point of view, the interest on central forces was related to 
the astronomical problem of planetary motions. However, there is no reason 
to consider them only under these circumstances. Another important issue 
that one can study within Classical Mechanics is the dispersion of particles 
by central forces. Of course, if the particles are of atomic size, we should keep 
in mind that the classical formalism may not give the right result because 
of quantum effects that begin to be important at those scales. Despite this, 
there are classical predictions that continue to be correct. Moreover, the 
main concepts of the dispersion phenomena are the same in both Classical 
Mechanics and Quantum Mechanics; thus, one can learn this scientific idiom 
in the classical picture, usually considered more convenient. 
In its one-body formulation, the dispersion problem refers to the action 
of the center of force on the trajectories of the coming particles. Let us 
consider a uniform beam of particles, (say electrons, protons, or planets and 
comets), but all of the same mass and energy impinging on a center of force. 
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We can assume that the force diminishes to zero at large distances. The 
incident beam is characterized by its intensity I (also called flux density), 
which is the number of particles that pass through per units of time and 
normal surface. When one particle comes closer and closer to the center of 
force will be attracted or repelled, and its orbit will deviate from the initial 
rectilinear path. Once it passed the center of force, the perturbative effects 
will diminish such that the orbit tends again to a streight line. In general, 
the final direction of the motion does not coincide with the incident one. 
One says that the particle has been dispersed. By definition, the differential 
cross section cr(0) is 

a{n)dn = (34) 

where dN is the number of particles dispersed per unit of time in the element 
of solid angle d^l around the O direction. In the case of central forces there 
is a high degree of symmetry around the incident beam axis. Therefore, the 
element of solid angle can be written 

dn = 2TTsmede, (35) 

where O is the angle between two incident dispersed directions, and is called 
the dispersion angle. 

For a given arbitrary particle the constants of the orbit and therefore the 
degree of dispersion are determined by its energy and angular momentum. 
It is convenient to express the latter in terms of a function of energy and 
the so-called impact parameter s, which by definition is the distance from 
the center of force to the straight suport line of the incident velocity. If uq 
is the incident velocity of the particle, we have 

I = muQS = s ■ ^/2mE. (36) 

Once E and s are fixed, the angle of dispersion is uniquely determined. 
For the time being, we suppose that different values of s cannot lead to 
the same dispersion angle. Therefore, the number of dispersed particles in 
the element of solid angle dVL between 9 and O + dQ should be equal to 
the number of incident particles whose impact parameter ranges within the 
corresponding s and s + ds: 

2ttIs \ds\ = 27rcr(e)/ sin e\de\. (37) 

In the equation (^^ we have introduced absolute values because while the 
number of particles is always positive, s and can vary in opposite direc- 
tions. If we consider s as a function of the energy and the corresponding 
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dispersion angle, 

s = s{e,E), 

the dependence of the cross section of Q will be given by 

.(6) ' '^^ 



sin B 



de 



(38) 



From the orbit equation (^), one can obtain directly a formal expression 
for the dispersion angle. In addition, for the sake of simplicity, we tackle 
the case of a pure repulsive dispersion. Since the orbit should be symmetric 
with respect to the direction of the periapsis, the dispersion angle is 



e = vr - 2^- 



(39) 



where ^ is the angle between the direction of the incident asymptote and 
the direction of the periapsis. In turn, ^ can be obtained from the equation 
( p5|) by making tq = oo when Oq = it (incident direction). Thus, 9 = tt — 
when r = the closest distance of the particle to the center of force. 
Then, one can easily obtain 



dr 



2/2mE 
1^ 



2mV _ 1 



Expressing / as a function of the impact parameter s (eq. 
is 

@ = TT-2 



sdr 



1 



E 



or 



e 



vr ■ 



sdu 



v{u) 



(40) 

the result 
(41) 

(42) 



The equations (^) and (42) are used rarely, as they do not enter in a direct 
way in the numerical calculation of the dispersion angle. However, when an 
analytic expression for the orbits is available, one can often get, merely by 
simple inspection, a relationship between and s. 



EXAMPLE : 

This example is very important from the historical point of view. It refers 
to the repulsive dispersion of charged particles in a Coulomb field. The field 
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is produced by a fixed charge —Ze and acts on incident particles of charge 
—Z'e\ therefore, the force can be written as follows 

that is, one deals with a repulsive inverse square force. The constant is 

k = -ZZe^. (43) 
The energy E is positive implying a hyperbolic orbit of eccentricity 



J 2EP 2 f 2Es X"^ 

where we have taken into account the equation (^) . If the angle 6' is taken 
to be vr, then from the equation (29) we come to the conclusion that the 
periapse corresponds to 6 = and the orbit equation reads 

1 mZZ'e^ r ^ -, . N 

- = ecos0-l. 45 

r 

The direction ^ of the incident asymptote is thus determined by the condi- 
tion r — > oo: ^ 

cos ^ = - , 
e 

that is, according to equation (^), 

. e 1 

sm — = — . 
2 e 

Then, 

2 2 
cot^ - = e2 _ 1^ 



and by means of equation (^) 



e lEs 
cot — 



2 ZZ'e^ ■ 

The functional relationship between the impact parameter and the disper- 
sion angle will be 

zz'e^ e 
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and by effecting the transformation required by the equation ( pq ) we find 
that a{@) is given by 

The equation (|47| ) gives the famous Rutherford scattering cross section de- 
rived by him for the dispersion of a particles on atomic nuclei. In the non- 
relativistic limit, the same result is provided by the quantum mechanical 
calculations. 

The concept of total cross section gt is very important in atomic physics. 
Its definition is 

aT= I (7{n)dn = 2tt r a{e)de . 

Jatt Jo 
However, if we calculate the total cross section for the Coulombian dispersion 



by substituting the equation (47) in the definition above we get an infinite 
result. The physical reason is easy to see. According to the definition, the 
total cross section is the number of particles per unit of incident intensity 
that are dispersed in all directions. The Coulombian field is an example of 
long-range force; its effects are still present in the infinite distance limit. The 
small deviation limit is valid only for particles of large impact parameter. 
Therefore, for an incident beam of infinite lateral extension all the particles 
will be dispersed and should be included in the total cross section. It is clear 
that the infinite value of cjt is not a special property of the Coulombian field 
and occurs for any type of long-range field. 



Further reading 

L.S. Brown, Forces giving no orbit precession, Am. J. Phys. 46, 930 (1978) 

H. Goldstein, More on the prehistory of the Laplace- Rung e-Lenz vector, Am. 
J. Phys. 44, 1123 (1976) 
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3. THE RIGID BODY 



Forward: Due to its particular features, the study of the motion of the 
rigid body has generated several interesting mathematical techniques and 
methods. In this chapter, we briefly present the basic rigid body concepts. 



CONTENTS: 

3.1 Definition 

3.2 Degrees of freedom 

3.3 Tensor of inertia (with example) 

3.4 Angular momentum 

3.5 Principal axes of inertia (with example) 

3.6 The theorem of parallel axes (with 2 examples) 

3.7 Dynamics of the rigid body (with example) 

3.8 Symmetrical top free of torques 

3.9 Euler angles 

3.10 Symmetrical top with a fixed point 
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3.1 Definition 

A rigid body (RB) is defined as a system of particles whose relative distances 
are forced to stay constant during the motion. 

3.2 Degrees of freedom 

In order to describe the general motion of a RB in the three-dimensional 
space one needs six variables, for example the three coordinates of the center 
of mass measured with respect to an inertial frame and three angles for 
labeling the orientation of the body in space (or of a fixed system within 
the body with the origin in the center of mass). In other words, in the 
three-dimensional space the RB can be described by at most six degrees of 
freedom. 

The number of degrees of freedom may be less when the rigid body is sub- 
jected to various conditions as follows: 

• If the RB rotates around a single axis there is only one degree of 
freedom (one angle). 

• If the RB moves in a plane, its motion can be described by five degrees 
of freedom (two coordinates and three angles). 

3.3 Tensor of inertia. 

We consider a body made of N particles of masses rua, a = 1,2,3..., N. 
If the body rotates at angular velocity a; around a fixed point in the body 
and this point, in turn, moves at velocity v with respect to a fixed inertial 
system, then the velocity of the ath particle w.r.t. the inertial system is 
given by 



Vet = V -I- a; X r, 



(1) 



The kinetic energy of the ath particle is 




(2) 



where 



Va = Va • Vq = (v (J X r^) • (v (J X r^) 
= V • V -I- 2v • (a; X Tq,) -|- (o; x Tq) • (a; x r^) 
= -v^ + 2v(a; x r^) + {cj x r^f. 



(4) 



(3) 
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Then the total energy is 

T = Ea^a = Ea I^T-aV^ +X;Q,mQ, [v- X Fq)] + 

If the origin is fixed to the sohd body, we can take it in the center of 
Thus, 

M 

and therefore we get 

T= ^Mv^ + i^m,, {uj X Taf 

a 



where 



1^ . ^2 



In Eq. ^ we use the vectorial identity 

(A X B)2 = A^B^ - (A • B)2 
to get the following form of the equation 

Trot = \^ma Uj'^r'^ - (t^ • 

which in terms of the components of uj and r 

(J = (a;i,a;2,t^3) and = (x^i, j:q,2, x^s) 
can be written as follows 

Now, we introduce 

= E '^ij'^i 

j 
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Trot = 2 miZ™'" [iOiUjdij (T,xi]}\ - UiUjXaiXaj\ 



Trot = 9 X] ^'■^j X] 



k 



We can write Trot as follows 



Trot — X^ lijl^i^j 



where 



k 



(10) 

(11) 

(12) 
(13) 



The nine quantities lij are the components of a new mathematical entity, de- 
noted by {lij} and called tensor of inertia. It can be written in a convenient 
way as a (3 X 3) matrix 

^ 111 Il2 Il3 ^ 

hi h2 hs = 
\ hi h2 h3 / 



{hj} = 

^ J2a "^a(a^a2 + ^as) 



^2a ^oiXalXa2 



\ 



- J2a maXalXa3 ^ 
J2a'^aXa3Xa2 J2a1T^a{xli + XI2) J 



X/Qf ^q(^q1 ~I~ X^ 



(14) 

We note that hj = Iji, and therefore is a symmetric tensor, implying 
that only six of the components arc independent. The diagonal elements of 
{hj} are called moments of inertia with respect to the axes of coordinates, 
whereas the negatives of the nondiagonal elements are called the products 
of inertia. For a continuous distribution of mass, of density p(r), {lij} is 
written in the following way 



•■13 



I P(r) 
Jv 



^ij X^ '^jt XiXj 
k 



dV. 



(15) 



EXAMPLE : 

Find the elements hj of the tensor of inertia {hj} for a cube of uniform 
density of side 6, mass M, with one corner placed at the origin. 

b b h 



'11 



P 



2 I 2 I 2 

X^ \ X'2 "r X^ 



XlXl 



dxidx2dx3 = p 



V 



+ x1)dxidx2dx. 



3 • 



000 
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The result of the three-dimensional integral is In = |(/96'^)^ = ^Mb"^. 

b b b 

h2 = J pi-xiX2)dV ^ ~P J J j {xiX2)dxidx2dx^ = -^pb^ = -^Mb^ . 

V 

We see that all the other integrals are equal, so that 



hi = I22 = h3 = ^Mb^ 



3 

J- 4 



leading to the following form of the matrix 



( ^Mb^ -iMb^ -iMb^ 
-^Mb'^ iMb'^ -{Mb"^ 
\ -\Mb'^ -\Mb'^ iMb'^ 



3.4 Angular Momentum 

The angular momentum of a RB made of particles of masses ma is given 

by 

L = ^r„ X p„ , (16) 

a 

where 

Pa = maVa = ma{uj X r^) . (17) 
Substituting ^ in (|l^), we get 



maYoL X X To 

a 



a 

Employing the vectorial identity 

A X (B X A) = (A • A)B - (A • B)A = A^B - (A • B)A 

leads to 

L =^mQ,(r^u; - Ya{ui ■ Va). 

a 

Considering the ith component of the vector L 
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and introducing the equation 



3 



we get 



Li = Ea (Ej ^ij^j Ek xlk) - (Ej 



— Ea Ej ^j^ij (Efc -^ak ^ai^aj) 
~ Ej Ea {^ij Ek ^^afc^^oi^^aj) 

Comparing with the equation (^) leads to 

Lj = ^ lijUJj . 
j 

This equation can also be written in the form 

L ={Iij}uj , 

or 





(18) 

(19) 
(20) 

(21) 

(22) 
(23) 



^ hi 1x2 hz 
hi h2 hs 
\ hi h2 hs 

The rotational kinetic energy, Trot, can be related to the angular momentum 
as follows: first, multiply the equation ( ^) by ^u>i 



-Li - -uJi^h 



(24) 



and next summing over all the i indices, gives 



ijLJiUJj . 



Compararing this equation with (|12D , we see that the second term is Trot- 
Therefore ^ ^ 

Trot = X! l^LiUJi = -L-uj . (25) 



Now, we substitute (22) in the equation (pq), getting the relationship be- 
tween Trot and the tensor of inertia 



Trot = -UJ- {hj} -UJ- 



(26) 
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3.5 Principal axes of inertia 

Taking the tensor of inertia {lij} diagonal, that is lij = Ii5ij, the rotational 
kinetic energy and the angular momentum are expressed as follows 

Trot = 2 ^ hj^i^j 

= - ^ 6ijIiUJiUJj 



Trot = -Y.I^^f (27) 

and 



2 



j 

L = luj. (28) 

To seek a diagonal form of {lij} is equivalent to finding a new system of 
three axes for which the kinetic energy and the angular momentum take the 
form given by ( p7|) and ([2^). In this case the axes are called principal axes 
of inertia. That means that given an inertial reference system within the 
body, we can pass from it to the principal axes by a particular orthogonal 
transformation, which is called transformation to the principal axes. 



Making equal the components of (22) and (P8|), we have 



Li = luji = hiuji + h2^2 + -^is^^s (29) 

L2 = IuJ2 = -^21^^! + 122^2 + hl,^^ (30) 
^3 = = hl^l + h2^2 + h3^3 ■ (31) 

This is a system of equations that can be rewritten as 

(111 - I)uJi + I12UJ2 + = (32) 

/21^1 + {I22 — 1)^2 + -^231^3 = 
IsiWl + /32W2 + (-^33 - -^)W3 = . 
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To get nontrivial solutions, the determinant of the system should be zero 



(111 - I)u)i /12W2 /l3^^3 
^21^1 {I22 — 1)^2 hs^-i 

I3ILOI h2<^2 (^33 — 1)1^3 



(33) 



This determinant leads to a polynomial of third order in I, known as the 
characteristic polynomial. The equation (^) is called the secular equation or 
characteristic equation. In practice, the principal moments of inertia, being 
the eigenvalues of I, are obtained as solutions of the secular equation. 

EXAMPLE : 

Determine the principal axes of inertia for the cube of the previous example. 

Substituting the values obtained in the previous example in the equation 
(H) we get: 



-3/3 (§/?-- 
-3/3 -3/3 



-3/5 
-3/3 



, 



where (3 = Mb^. Thus, the characteristic equation will be 



12'^ 



The solutions, i.e., the principal moments of inertia are: 



^1 = 





3 -^2'^ 



whose corresponding eigenvalues are given by 
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The diagonalized {lij} will be 







\ 






12 



3.6 The theorem of parallel axes 

We suppose that the system xi,X2,X3 has the origin in the center of mass 
of the RB. A second system Xi, X2, X^, has the origin in another position 
w.r.t. the first system. The only imposed condition on them is to be parallel. 
We define the vectors r = {xi,X2,X3), R = {Xi, X2, X^) y a = (01,02,03) 
in such a way that R = r + a, or in component form 



Xi — Xi -\- di. 



(34) 



Let Jij be the components of the tensor of inertia w.r.t. the system X1X2X3, 



rrir 



k 



(35) 



We substitute (B3) in (pg), 



J, 



E 



k 



= [Ea "T-a i^ij Y^ki^^kf " XaiXaj)] + J2a {Sij J2k 4 " O-ittj) (36) 
+ Efc 2afc(5ij (X;^ rriaXak) - aj {J2a f^aXaj) - ai {J2a ITT-aXai)] ■ 

Since the center of mass coordinate is defined as 

we take into account that we have already set the origin in the center of 
mass, i.e., 

(xi,X2,X3) = (0,0,0) . 

Now, if we also compare the first term in (^) with the equation ([l^), we 
have 

Jij = lij + M{a^6ij - ttiaj) (37) 
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and therefore the elements of the tensor of inertia lij for the center of mass 
system will be given by: 



Jij — M{5ija^ — aiUj) . 



(38) 



This is known as the theorem of the parallel axes. 



EXAMPLE : 

Find lij for the previous cube w.r.t. a reference system parallel to the system 
in the first example and with the origin in the center of mass. 
We already know from the previous example that: 



i/3 



1/3 -IP 



Now, since the vector a = (2; 2' 2. 
(H) and the fact that (3 = Ml? to get, 



and a? = we can use the equation 



Therefore 



hi 
h2 

h2 
h2 



{1} 



a?) 
4) 



Jii - M(a2 

J22 - M(a2 - al) = \Mb'^ 
J33 - M(a2 - al) 



1Mb'' 



J12 - M{-aia2) = 
= -^13 = -^23 = . 



6 



V 




(39) 
(40) 
(41) 
(42) 
(43) 



EXAMPLE . 

We consider the case for which the vector a 
the new tensor of inertia will be: 



hi 
I22 
hs 
I12 



J 



11 



M(a2 



(0: 



M 



b b- 
2' 2' 



62 



and = ^. Then, 



J22 - M(a2 - al) = [1Mb') - M - f 
J33 _ Mia' - al) = {1Mb') - M - f 
= J12 - M{-aia2) = {-jMb') - M(0) = 



= M' 

D 

= ^Mb' 
= l^Mb' 
-1Mb' 



(44) 
(45) 
(46) 
(47) 
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Ii3 = Ji3-M(-aia3) = (-iM62 
/23 = ^23 - M(-a2a3) = (-iM62 

It follows that {lij} is equal to: 



( iMh'' 
-\M}? 
-\M}? 



-\Mb^ 




M(0) = - 
-M{\Mb'^ 



iMb^ 
1 = 0. 



-\Mb^ 




(48) 
(49) 



3.7 The dynamics of the rigid body 

The rate of change in time of the angular momentum L is given by: 



= ^^^\ (50) 

inertial 



For the description w.r.t. the system fixed to the body we have to use the 
operator identity 



Applying this operator to the equation (|5 



+ UJX . (51) 

body 



+ CJXL. (52) 

body 



Then, instead of (^) we shall have 

+ wxL = N. (53) 



dt 



body 



Now we project the equation (^) onto the principal axes of inertia, that we 
call {xi,X2,xz); then Tj-ot and L take by far simpler forms, e.g., 

Li = liLVi . (54) 

The ith component of (]5^ ) is 

^ + eijkUJjLk = Ni . (55) 



42 



Now projecting onto the principal axes of inertia and using the equation 
(|5^, one can put (|55|) in the form: 

du>i 

li — + eijkUJjUJkh = Ni (56) 

since the principal elements of inertia are independent of time. Thus, we 
obtain the following system of equations known as Euler's equations 

IliOl + UJ2^^3{l2 - h) = ^1 (57) 
I2UJ2 + {I3 - h) = N2 
I3LU3 + uJiuJ2{h - h) = N3 . 



EXAMPLE : 

For the rolling and sliding of a billiard ball, prove that after a horizontal 
kick the ball slips a distance 

X, - 



where uq is the initial velocity. Then, it starts rolling without gliding at the 
time 

= 

SOLUTION: When the impulsive force stops, the initial conditions are 

XO = 0, ±0 = Uq 

(f) = 0, (/) = . 

The friction force is given by 

F/ = -figei , 
and the equation of motion reads 

X = -figM. (58) 

The equation for L is 

^=I,4> = N3 (59) 
at 



43 



where I-\ is 



and 



P(r) 



2 2 
Xi — X2 



dx\dx2dx2 = —Ma 
5 



Substituting I3 and A^s in (p^), one gets 



(60) 



Integrating once both (58) and (|60|), we get 

i; = -/i^ft + Ci (61) 

ac/) = ^igt + C72 • (62) 
To these equations we apply the initial conditions to put them in the form 

x(t) = —fj.gt + no (63) 

(64) 

The condition of pure rolling (no friction) is 

x{t) = a^{t). (65) 
From (^) and ( |65|) evaluated at ti we get 



a0(i) = 2 W*- 



-ugti = -figh + uo 



2uo 



(66) 



Now we integrate ( p3[ ) once again and applying the initial conditions we get 



x(t) = -pt^^ + ^^0* • 
By evaluating (^) and (^3|) at time ti we are led to 



(67) 
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3.8 Symmetrical top free of torques 

A symmetric top is any solid of revolution. If the moments of inertia are 

h = h = I3 spherical top 

h = I2 7^ I3 symmetric top 

h 7^ h 7^ I3 asymmetric top. 

Let us consider the symmetric top Ii = I2 ^ I^- In this case the axis 
is the axis of symmetry. The Euler equations projected onto the principal 
axes of inertia read 

Iiuji + oj20J3{h - h) = Ni (68) 
l2io2 + 0J3UJi{h - h) = N2 (69) 
+ uJiUJ2{Ii - h) = N3. (70) 
Since the system we consider here is free of torques 

Ni = N2 = Ns = , (71) 



we use Ii = I2 in (71) to get 

/it^i+a;2a;3(/2-/3) = (72) 

I2UJ2 + uj3u;i{l3 - /i) = (73) 

hdJ3 = 0. (74) 



The equation (74) implies that 

UJ3 = const. 



The equations (72) and (^) are rewritten as follows: 

uji = —nuj2 where Q = UJ3 ( — ) (75) 



UJ2 = — Ocji . (76) 
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Multiplying (|7^) by i and summing it to (|75|), we have 

{uJi+iL02) = i^{uJi + i0J2)- 
If we write r]{t) = uJi{t) + iuj2{t), then 

r](t) - inr]{t) = . 

The solution is 



This implies 
Thus, 



7/(t) = Aexp{ii}t) . 
{loi + iuj2) = A cos(ilt) + i sin(Ot) 



uJi = Acos{nt) (77) 
0J2 = Asm{nt). (78) 
The modulus of the vector uj does not change in time 



LO = UJ 



\Jlo\ + a;2 + W3 = \7 + Wo = consi 



This vector performs a precessional motion of precession frequency given 
by 

Moreover, we notice that is constant. 

If we denote by A the angle between uj and the equations ([77[) and (|7: 
take the form 

wi = cj sin A cos(f]t) 
L02 = sin A sin(r2f) 
wi = a; cos A , 

where A = lo sin A. 

For a flattened body of revolution Ii = I2 = I12 and I3 > Ii. For example, 
for the case of the earth 

n fh-h2\ W3 
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The observations point to a mean value of fourteen months ~ 450 days. 
(This is due to the fact that the earth is not strictly a RB; there is also an 
internal liquid structure). 

3.9 Euler angles 

As we already know, a rotation can be described by a rotation matrix A by 
means of the equation 

X = Ax' . (79) 

X represents the set of axes of the system rotated w.r.t. the system whose 
axes are represented by x'. The rotation A can be accomplished through 
a set of "partial" rotations A = AiA2...An. There are many possibilities to 
choose these A's. One of them is the set of angles (j),6 and ip called Euler 
angles. The partial rotations are in this case the following: 

• A rotation around the X'3 axis of angle f (in the positive trigonometric 
sense). The corresponding matrix is: 



A(^ 



A rotation of angle 9 around the X'l axis (positive sense). The associ- 
ated matrix is: 

'10 

Xg = \ cos 9 sin ( 
— sin 9 cos ( 

• A rotation of angle (p around the X3 axis (positive sense) ; the assoiated 
matrix is: 

(coscj) sin</) ^ 
— sin (j) cos (j) 
01; 

The full transformation of the system of axes {X'l, X'2, X'^} to the system 
of axes {Xi, X2, X^} is given by (^), where 

Doing the product of matrices, we get 

All = cos if cos (p — cos 9 sin (p sin ip 
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where 



A21 = — sin {p cos (j) — cos sin (j) cos 
A31 = sin 9 sin 
A12 = cos sin (j) + cos cos cf) sin 
A22 = — sin ip sin (f) + cos 6* cos (j) sin 
A32 = — sin cos 
A13 = sin ip cos (j) 
A23 = cos sin 9 
A33 = cos 9 

All A12 Ai3 
A = I A21 A22 A23 
A31 A32 A33 



Now, we take into account that: 

• is along the (fixed) axis. 

• ^ is along the so-called line of nodes. 

• is along the X3 axis (of the body). 



This allows to write the three components of each of the three vectors in the 
system {Xi,X2,X^^ as follows: 

4>i = (j)s\n9 svmp^ 9i = 9cosip cpi = 

(p2 = <psin9 cos ip, 92 = — 9 simp (^2 = 

4>3 = ^cos9, 6*3 = (p3 = (p . 

Then, 



LO = (p + 9 + ip 

= (J>1 + ^1) , (J>2 + 92 + <P2) , (J>3 + O3 + 933) 

Thus, we are led to the following components of u: 

uj\ = ^ sin 9 simp + 9 cos (p 
U2 = 4> sin 9 cos ip — 9 sin ip 
0^3 = cos 9 + (p. 
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3.10 Symmetrical top with a fixed point 

As a more complicated example of the methods used to describe the dynam- 
ics of the rigid body, we shall consider the motion of a symmetric body in a 
unifom gravitational field when a point of the axis of symmetry is fixed in 
the space. 

The axis of symmetry is of course one of the principal axes, and we shall 
take it as the z axis of the body-fixed reference system. Since there is a 
fixed point, the configuration of the top will be determined by the three 
Euler angles: Q measuring the deviation of z from the vertical, ^, giving the 
azimuth of the top w.r.t. the vertical, and which is the rotation angle of 
the top w.r.t. its proper z. The distance from the center of gravity to the 
fixed point will be denoted by I. To get a solution to the motion of the top 
we shall use the method of Lagrange instead of the Euler equations. 
The kinetic energy is: 

or, in terms of the Euler angles: 

T = gij^2 0^02^ + ^j^(^^cose + iff . 

According to an elementary theorem, in a constant gravitational field the 
potential energy of a body is the same with that of a material point of equal 
mass concentrated in its center of mass. A formal proof is as follows. The 
potential energy of the body is the sum of the potential energies of all its 
particles: 

V = -ruiVi ■ g , (80) 

where g is the constant acceleration of gravity. According to the definition 
of the center of mass, this is equivalent to 

V = -MR, • g, (81) 

thus proving the theorem. The potential energy is a function of the Euler 
angles: 

V = Mgl cose, (82) 

and the Lagrangian will be 

L = ^hi^"^ sin^ e + e'^) + ^hi^ cos e + tpf -Mgl cose. (83) 
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We note that (p and (p are cyclic coordinates, and therefore and p^p are 
constants of motion. 

dL 

Pip = -qT = Isi^p + (pcos 0) = const (84) 

and 

dL 

ps = — - = Ii(j)s\v? 9 + I^{cf) cos^ 9 + if) cos 9) = const. (85) 



From the equation (p3) we get 



that we substitute in m 



Pip -h<i> cos 9 

LP = , (86) 

h 



T X ■ a ^ T (X , Pv - h4>cos9 
p^ = — r=ii0sm 9 + l3[(p COS 9 -\ — cos9)= const. 

p^ = Ii4>siv? 9 + p^cos9 , 

where from we get 

cosg_ ^^^^ 
Ji sin 9 



Substituting it in (86) one gets 



Pip P<t> - Pip cos ( 

(p = -r- ^-75 cos ( 

^ Is h sin2 9 



Now, since the system is conservative, another integral of motion is the 
energy 

E = T + V = hi{^'^ sin^ 9 + 9"^) + h^i^ cos 9 + pf + Mgl cos 9. 

The quantity I3UJ3 = p^ is an integral of motion. Multiplying this constant 
by p^uj3 we get 

hPipi^s = pI^^s 
li^i = p1^3 

2 ^ 2 I3 
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The quantity 5/3^3 is a constant. Therefore, we can define the quantity 



E' 



1 2 
E — -huJ^ = const. 

^hO^ + sin^ 6 + Mgl cos 6 , 



wherefrom we can identify 



V{e) = -<j)'^sm^e + Mgl cost 



V{9) = -h 



Pip cos 9 



Thus, ^is: 



h sin2 9 
1 



sin2 6' + M5Z cos 6*. 



K= 11^9^ + V{9) . 



From this equation we get 9 = 



1/2 



(89) 



{E- V{9)) , which leads to 
d9 



(^) {E-V{9)) 



(90) 



Performing the integral in (|90D one gets t = f{9), and therefore, in principle, 
one can get 9{t). Then, 9{t) is replaced by (p and ip (in eqs. ( |87| ) and (|88|)) 
and integrating them we can obtain the complete solution of the problem. 
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4. SMALL OSCILLATIONS 



Forward: A familiar type of motion in mechanical and many other systems 
are the small oscillations (vibrations). They can be met as atomic and 
molecular vibrations, electric circuits, acoustics, and so on. In general, any 
motion in the neighborhood of stable equilibria is vibrational. 



CONTENTS: 

4.1 THE SIMPLE HARMONIC OSCILLATOR 

4.2 FORCED HARMONIC OSCILLATOR 

4.3 DAMPED HARMONIC OSCILLATORS 

4.4 NORMAL MODES 

4.5 PARAMETRIC RESONANCE 
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4.1 THE SIMPLE HARMONIC OSCILLATOR 



A system is at stable equilibrium when its potential energy U (q) is at min- 
imum; when the system is slightly displaced from the equilibrium position, 
a force —dU/dq occurs which acts to restore the equilibrium. Let q^ be the 
value of the generalized coordinate corresponding to the equilibrium posi- 
tion. Expanding U (q) — U{qo) in a Taylor series of (7 — go for small deviations 
from the equilibrium 



This means that there are no external forces acting on the system and the 
zero has been chosen at the equilibrium position; moreover, higher-order 
terms have been neglected. The coefficient k represents the value of the 
second derivative of U{q) for q=qo- For simplicity reasons we denote 

x = q-qo 

for which the potential energy can be written as: 




donde: 



dU 
dq 
U{q) 



= 



= 





For simplicity reasons we denote 



X = q 



for which the potential energy can be 



written as: 





The kinetic energy of a system is 




and using (|2[) and we get the Lagrangian of a system performing linear 
oscillations (so-called linear oscillator): 
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1 -2 1 2 

L = -m X kx . (4) 

The equation of motion corresponding to this L is: 

m X +kx = , 

or 

X +w'^x = , (5) 

where w'^ = \fkjra. This differential equation has two independent solu- 
tions: cos wt and sinwt, from which one can form the general solution: 

x = ci cos wt + C2 sin wt , (6) 
or, we can also write the solution in the form: 

X = acos{wt + a) . (7) 

Since cos{wt + a) = cos wt cos a — sinwtsma, by comparing with (^), one can 
see that the arbitrary constants a and a are related to ci and C2 as follows: 



(cf+cl), y tana = — ci/c2 . 

Thus, any system in the neighborhood of the stable equilibrium position 
performs harmonic oscillatory motion. The a coefficient in (^) is the am- 
plitude of the oscillations, whereas the argument of the cosine function is 
the phase of the harmonic oscillation; a is the initial value of the phase, 
which depends on the chosen origin of time. The quantity w is the angular 
frequency of the oscillations, which does not depend on the initial conditions 
of the system, being a proper characteristic of the harmonic oscillations. 
Quite often the solution is expressed as the real part of a complex quantity 

X = Re [Aexp{iwt)] 

where A is the complex amplitude, whose modulus gives the ordinary am- 
plitude: 

A = oexp(za) . 
The energy of a system in small oscillatory motion is: 
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or by substituting (|^ 



^ 1 -2 1, 2 

E = -m X H — kx , 
2 2 



1 2 2 

£/ = —mw a . 



Now, we consider the case of n degrees of freedom. In this case, taking the 
sum of exterior forces as zero, the generahzed force will be given by 

dU 

Qi = -^ = 0. (8) 
dqi 

Repeating the procedure for the case of a single degree of freedom, we expand 
the potential energy in Taylor series taking the minimum of the potential 
energy at qi = qiQ. Introducing small oscillation coordinates 

Xi = qi- qio , 
we can write the series as follows 



Uiquq2, qn) = U{q,o, <?20, gno)+E (1^) ^ ^^ + ^ E ^ ^^^. + - 

(9) 

Under the same considerations as given for (^), we obtain: 

U{qi,q2, q-n) = U =^Y1 hjXiXj . (10) 

From (^) one notes that kij = kji, i.e., they are symmetric w.r.t. their 
subindices. 

Let us look now to the kinetic energy, which, in general, is of the form 

— ajj(5) XiXj , 

where the functions of the coordinates only. Denoting them by 

ttij = niij the kinetic energy will be 

T= \ ^'mij XiXj . (11) 
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We can pass now to the Lagrangian for the system of n degrees of freedom 



L — T U — 2 ^ y (^l7 Xj^Xj kijXiXj^ 



(12) 



This Lagrangian leads to the following set of simultaneous differential equa- 
tions of motion 



d dL dL 



dxi 



= 



or 



didxi 
'^{ruij Xj +kijXj) = . 



(13) 



(14) 



This is a linear system of homogeneous equations, which can be considered 
as the n components of the matricial equation 



(M)(X) + (K)(X)=0, 
where the matrices are defined by: 



(M) 



/ mil "112 

17121 17122 



min \ 
m2n 



{K) 



\ rn-ai 


mn2 ■ 


TTlnn 


( kn 


ki2 ■■ 


■ kin ^ 


k2i 


k22 ■■ 


■ k2n 



\ kji\ kn2 •■• knn / 
( XI \ 



X2 



{X) 



/ XI \ 

X2 

\ Xn J 



(15) 



(16) 



(17) 



(18) 



(19) 
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Similarly to the one dof system, we look for n unknown functions Xj{t) of 
the form 



Ajex.p{iwt) , (20) 



where Aj are constants to be determined. Substituting (pO|) in (|14D and 
dividing by exp{iwt), one gets a linear system of algebraic and homogeneous 
equations, which should be fulfilled by Aj. 

Y.i-w'^mik + kik)Ak = . (21) 
j 

This system has nponzero solutions if the determinant of its coefficients is 
zero. 



^ = 0. (22) 



This is the characteristic equation of order n w.r.t. w'^. In general, it has n 
different real and positive solutions Wa (a = 1,2, ...,n). The Wa are called 
proper frequencies of the system. Multiplying by A* and summing over i 
one gets 

^{-w'^rriij + kij)A*Aj = , 



where from 



y ] kijA^ Ail ^ ^ niijA^ Ai . 



Since the coefficients kij and niij are real and symmetric, the quadratic forms 
of the numerator and denominator are real, and being essentially positive 
one concludes that w'^ are equally positive. 
EXAMPLE 

As an example we model the equations of motion of a double pendulum. 
The potential energy of this system with two degrees of freedom is 

U = migli{l — cos 6*1) + m2gli{l — cos^i) + m2gl2{^ — cos 6*2) . 
Applying (^, one gets 

U = ^(mi + m2)gh6l + ]^m2gl20l . 
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Comparing with (p!o|), we identify 



hi = {'mi + fn2)li 
kl2 = ^21 = 

k22 = rn2gl2 ■ 



For the kinetic energy one gets 



1 -21 -2 

T = -{nil + m2)ll Oi +i^rn2ll 02 +"1-2^1^2 6*1612 ■ 



Identifying terms from the comparison with (|ll|) we find 

mil = {mi+m2)ll 
17112 = 17121 = m-2^1^2 
m22 = m2ll. 

Substituting the energies in (|l2|) one obtains the Lagrangian for the double 
pendulum oscillator and as the final result the equations of motion for this 
case: 

I mil mi2 \ f 61 \ ( ku \ / ^1 \ g 
1^ m2i ^22 y \^ ^2 / V k22 j \02 j 

4.2 FORCED HARMONIC OSCILLATOR 

If an external weak force acts on an oscillator system the oscillations of the 
system are known as forced oscillations. 

Besides its proper potential energy the system gets a supplementary po- 
tential energy Ue{x,t) due to the external field. Expanding the latter in a 
Taylor series of small amplitudes x we get: 



Ue{x,t) ^ Ue{0,t)+X 



dUe 

dx 



x=0 



The second term is the external force acting on the system at its equilibrium 
position, that we denote by F{t). Then, the Lagrangian reads 

1 .2 1 „ 
L = -mx --kx^ + xF{t). (23) 

The corresponding equation of motion is 
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m X +kx = F{t) , 



or 



X +w'^x = F{t)/m , (24) 

where w is the frequency of the proper oscillations. The general solution of 
this equation is the sum of the solution of the homogeneous equation and a 
particular solution of the nonhomogeneous equation 

X — Xfi -\- Xp . 

We shall study the case in which the external force is periodic in time of 
frequency 7 forma 



= /cos(7t + /3) 



The particular solution of (|2J) is sought in the form xi = bcos{'yt + (5) and 
by substituting it one finds that the relationship h = f /m{vP' — 7^) should 
be fulfilled. Adding up both solutions, one gets the general solution 

X = a cos{wt + a) + f /■m{w'^ — ^'^) cos{'jt + /S) . (25) 

This result is a sum of two oscillations: one due to the proper frequency and 
another at the frequency of the external force. 

The equation (^) can in general be integrated for an arbitrary external 
force. Writing it in the form 

— {x +iwx) — iw{x +iwx) = — F{t) , 
at m 

and making ^ =x +iwx, we have 

— ^ — iw^ = F{t)/m . 
The solution to this equation is ^ = A(t) exp{iwt); for A(t) one gets 

A= F{t)ex];){—iwt)/m . 
Integrating it leads to the solution 
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/■* 1 

^ = cxpUwt) / —F(t) exp(-iwt)dt + Co ■ (26) 
JO in 

This is the general solution wc look for; the function x(t) is given by the 
imaginary part of the general solution divided by w. 
EXAMPLE 

We give here an example of employing the previous equation. 

Determine the final amplitude of oscillations of a system acted by an cxtcnal 

force Fq = const, during a limited time T. For this time interval we have 

£ = — cyipUwt) / ex.p(—iwt)dt , 
m Jo 

f = [1 — exp(—iwt)] expUwt) . 

iwm 



Using 1^1 = a^w"^ we obtain 



a = ^ sm( -'u;T) 

mw^ 2 



4.3 DAMPED HARMONIC OSCILLATOR 

Until now we have studied oscillatory motions in free space (the vacuum), 
or when the effects of the medium through which the oscillator moves are 
negligeable. However, when a system moves through a medium its motion 
is retarded by the reaction of the latter. There is a dissipation of the energy 
in heat or other forms of energy. We are interested in a simple description 
of the dissipation phenomena. 

The reaction of the medium can be imagined in terms of friction forces. 
When they are small we can expand them in powers of the velocity. The 
zero-order term is zero because there is no friction force acting on a body 
at rest. Thus, the lowest order nonzero term is proportional to the velocity, 
and moreover we shall neglect all higher-order terms 

fr~ 2; , 

where x is the generalized coordinate and a is a positive coefficient; the 
minus sign shows the opositc direction to that of the moving system. Adding 
this force to the equation of motion we get 

m x= —kx — ax , 
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or 



x= —kx/m — ax /m . (27) 

Writing k/m = and a/m = 2A; where Wo is the frequency of free oscilla- 
tions of the system and A is the damping coefficient. Therefore 

x +2A x +wlx = . 

The solution of this equation is sought of the type x = exp(rt), which we 
substitute back in the equation to get the characteristic equation for r. Thus 

+ 2A + u;^ = , 

where from 

We are thus led to the following general solution of the equation of motion 

X = ci exp(rii) + C2 exp(r2i) . 

Among the roots r we shall look at the following particular cases: 
(i) A < Wo- One gets complex conjugate solutions. The solution is 

X = Re l^exp 



Xt + iJiw^o-X^) 



where A is an arbitrary complex constant. The solution can be written of 
the form 



= aexp(— At) cos{wt + a), where w = ^ [Wg — A^) , 



(28) 



where a and a are real constants. Thus, one can say that a damped oscil- 
lation is a harmonic oscillation with an exponentially decreasing amplitude. 
The rate of decreasing of the amplitude is determined by the exponent A. 
Moreover, the frequency w is smaller than that of free oscillations. 

(ii) A > Wo- Then, both r are real and negative. The general form of the 
solution is: 



X = ci exp 



A - ^(A2 - w^J 



t> + C2 exp 



t 
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If the friction is large, the motion is just a monotone decaying amplitude 
asymptotically (t oo) tending to the equilibrium position (without oscil- 
lations). This type of motion is called aperiodic. 

(iii) X = Wo- Then r = — A, whose general solution is 

X = (ci + C2t) exp(— At) . 

If we generalize to systems of n degrees of freedom, the generalized fric- 
tion forces corresponding to the coordinates Xi are linear functions of the 
velocities 



fr,i — ^ ^ Oiij ■ (^9) 
j 

Using Oik = aki, one can also write 

dF 

Ir,i ri • ) 

d Xi 

where F = cnj ^i^j is called the dissipative function. The differential 

equation is obtained by adding up all these forces to ([l^ 

^(mjj Xj +kijXj) = - ^ aij Xi . (30) 

j 

Employing 



Xk = Ak exp(rt) 

in ( ^0| ) and deviding by exp(rt), one can obtain the following system of linear 
algebraic equations for the constants Aj 

'^{rriijr^ + aijr + kij)Aj = . 



Making equal to zero the determinant of this system, one gets the corre- 
sponding characteristic equation 



rriijr + Oijr + kij 
This is an equation for r of degree 2n. 



. 



(31) 
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4.4 NORMAL MODES 



Before defining the normal modes, we rewrite (15) as follows 



M X) + K\X) 







where \X) is the n-dimensional vector whose matrix representation is (|TW); 
M and K are two operators having the matrix representation given by (|lq ) 
and (p!7|), respectively. We have thus an operatorial equation. Since M is 
a nonsingular and symmetric operator, the inverse operator M^^ and the 
operators M"^/^ and M~^/^ are well defined. In this case, we can express the 
operatorial equation in the form 



dt2 



M^/^ \x) = -M-^/^KM-^/^M^/^ \X) 



or more compactly 



where 



X 



X 



(32) 



x) = Ml/2 ij^v 



and 



Since M"^/^ ^j^d K are symmetric operators, then A is also symmetric. If 
we use orthogonal eigenvectors as a vectorial base (for example, the three- 
dimensional Euclidean space) , the matrix representation of the operator can 
be diagonal, e.g.. 



Xi6ij 



Let us consider the following eigenvalue problem 

^ \Pi) = \Pi) ) 

where \pi) is an orthogonal set of eigenvectors. Or 



(33) 
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The eigenvalues are obtained by multiplying both sides by leading to 

{pi\M-y^KM-^/^\p,) 



A,: 



{Pi \Pi 



Since the potential and kinetic energies are considered positive quantities, 
one should take 



(pi|M-i/2KM-i/2|p.))o 



and therefore 



This leads to the set 



A ,; > . 



Ai = wf . 



If we express the vector in terms of these eigenvectors of A, 

x)=Y.yi\x 

i 

where 

Vi = {pi X 

Inserting this result in the equation of motion ( |32D , we obtain 



(34) 



dt2 ^ 



The scalar product of this equation with the constant eigenvector {pj\ leads 
to the equation of motion for the generalized coordinate yj 



The solution of this equation reads 

yj = Aj cos{wjt + 



(35) 
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Use of these new generalized harmonic coordinates lead to a set of indepen- 
dent equations of motion. The relationship between yj and Xj is given by 

r 



Vj = Pjl ^1 +/'i2 X2 +... + Pjn Xn ■ 

The components pji {I = 1, 2, .., n) are determined by solving the eigenvalue 
problem given by (^). The new coordinates are called normal coordinates 
and the Wj are known as the normal frequencies. The equivalent matrix 
form (^5|) is 



™(i) 

^2 



Aj cos^Wjt + 



/ Pji \ 

Pj2 
V PJn J 



(36) 



These are the normal vibrational modes of the system. One reason for 
introducing the coordinates yj is found from the expression for the kinetic 
energy, which is seen to be invariant under the rotation to the new axes. 

1 



EXAMPLE 



Apply the matricial procedure as already shown, given the following equa- 
tions of motion 





Comparing with (^), we identify the operator A. To find the eigenvectores, 
we use (^3|) getting 




'3 / 

The characteristic equation for Aj is 

det(A - Xil) = , 
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and by substituting the values 



5 - A 1 

2- A 

1 5- A 



Solving the equation one gets Aj = 2, 4, 6. For A = 4 

( Pi\ 

we have the following set of equations 



5 1 

2 

1 5 



Pi 

P2 
P3 



(5-4)pi+/93 = 
2p2 - 4p2 = 
Pi + (5 - 4)p3 = . 

Taking into account the normalization condition, one is led to the following 
values 



Pi 

P2 



Therefore 



\p\= 



By the same means one gets 



-P3 



V2 



. 



|PA=6) 
|PA=2) 



1 

V2 


1 





/ 1 



V 1 
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Thus, the new vectorial space is determined by 



^ ^ 

1 
_J_ J_ 



where from 



{pi 



Thus, the normal coordinates are given by (|34 






4.5 PARAMETRIC RESONANCE 

The important phenomenon of parametric resonance shows up for systems 
initially at rest in un unstable equilibrium position, say x = 0; thus, the 
slightest deviation from this position produces a displacement growing rapidly 
(exponentially) in time. This is different from the ordinary resonances, 
where the displacement grows only linearly in time. 

The parameters of a linear system are the coefficients m and k of the La- 
grangian (Q); if they are functions of time, the equation of motion is: 

^(mx) + A;x = 0. (37) 
If we take a constant mass, the previous equation can be written in the form 

-^+w^{t)x = 0. (38) 

The function w{t) is given by the problem at hand. Assuming it a periodic 
function of frequency 7 (of period T = 27r/7), i.e., 

w{t + T) =w{t) , 

any equation of the type ( ^8|) is invariant w.r.t. the transformation t t+T. 
Thus, if x{t) is one of its solutions, x{t + T) is also a solution. Let xi{t) and 
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X2{t) be two independent solutions of 4.6.2). They should change to itselves 
in a linear combination when t ^ t + T. Thus, one gets 

xi{t + T) = fiix{t) (39) 

X2{t + T) = fl2xit) , 

or, in general 

x^{t) = t/fF{t) 

X2{t) = /x^Glt) , 

where F{t) and G{t) are periodical functions in time of period T. The 
relationship between these constants can be obtained by manipulating the 
following equations 

xi +w'^{t)xi = 
X2 +w^{t)x2 = . 

Multiplying by X2 and xi, respectively, and substracting term by term, we 
get 

d . 

Xi X2- X2 Xi = —iXi X2- X2 Xi) = , 

at 

or 

x'l X2— X2 xi = const. . 

Substituting t by t + T in the previous equation, the right hand side is 

multiplied by /xi^2 (see eqs. (|39|)); thus, it is obvious that the following 
condition holds 

/xi/i2 = 1 , (40) 

where one should take into account (|38| ) and the fact that the coefficients 
are real. If x{t) is one integral of this equation, then x*{t) is also a solution. 
Therefore /xi, /X2 should coincide with fil, This leads to either fii=^i2 
/ii and ^2 both real. In the first case, based on ( ^0[ ) one gets /ii = 1/ //|, 
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1 2 

that is equivalent to = |/i2 
are of the form 



2 



1. In the second case, the two solutions 



X2{t) 



li-'/^G{t) . 



One of these functions grows exponentially in time, which is the character- 
istic feature of the parametric resonance. 



* H. Goldstein, Classical mechanics, Second ed. (Addison- Wesley, 1981). 

* L. D. Landau k. E. M. Lifshitz, Mechanics, (Pergammon, 1976). 

* W. Hauser, Introduction to the principles of mechanics, (Wesley, 1965). 

* E.I. Butikov, Parametric Resonance, Computing in Science &; 
Engineering, May/June 1999, pp. 76-83 (http://computer.org). 
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5. CANONICAL TRANSFORMATIONS 



Forward: The main idea of canonical transformations is to find all those 
coordinate systems in the phase space for which the form of the Hamilton 
eqs is invariant for whatever Hamiltonian. In applications one chooses the 
coordinate system that allows a simple solution of the problem at hand. 

CONTENTS: 

5.1 Definitions, Hamiltonians and Kamiltonians 

5.2 Necessary and sufficient conditions for canonicity 

5.3 Example of application of a canonical transformation 
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5.1 Definitions, Hamiltonians and Kamiltonians 



For the time-independent and time-dependent cases, respectively, one de- 
fines a canonical transformation as follows 

Definition 1: A time-independent transformation Q = Q{q,p), and P = 
P{q,p) is called canonical if and only if there is a function F{q,p) such that 

dF{q,p) = ^Pidqi -^Pi{q,p)dQi{q,p) . 

i i 

Definition 2: A time-dependent transformation Q = Q{q,p,t), and P = 
P{q,p,t) is called canonical if and only if there is a function F{q,p,t) such 
that for an arbitrary fixed time t = to 

dF{p,q,to) = ^Pidqi - ^Pi{q,p,to)dQi{p,qM) , 
i i 

where 

^TPf +\ sr^dF{p,q,to) ^yr^ dF{p,q,to) 
dF{p, q, to) = 2^ — dqi + 2^ — dpi 

and 

+^ ST- dQ{p,q,to) , to) , 
dQ{p, q, to) = 2^ dqi + dpi 

Example: Prove that the following transformation is canonical 

Q = Tan'^ . 

Solution: According to the first definition we have to check that pdq — PdQ 
is an exact differential. Substituting P and Q in the definition we get 

^ Dj/n J ^f2. 2^pdq-qdq , f pq\ 

pdq-PdQ=pdq--{p +q) ^2^q2 = [y ) ' 

We can see that indeed the given transformation is canonical. We know 
that a dynamical system is usually characterized by a Hamiltonian H = 
H{q,p,t), where q = q{qi,q2, ■■■,qn), and p = p{pi,P2, ■■■,Pn)- Therefore, 
the dynamics of the system fulfills a set of 2n first-order differential eqs 
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(Hamilton's eqs.) 

dH 

dpi 
dH 

Let us denote the coordinate transformations in the phase space by 

Qj=Qj{q,P,t) (3) 
Pj = Pj{q,p,t) . (4) 

According to the aforementioned principle for the set of canonical trans- 
formations denoted by (3) and (4), analogously to (1) and (2), there is a 
function K = K{Q, P, t) such that we can write 

The relationship between the Hamiltonian H and the Kamiltonian can 
be obtained arguing as follow^. 

According to Hamilton's principle, the real trajectory of a classical system 
can be obtained from the variation of the action integral 

6 JiY.p^dq,-Hdt)=0 . (7) 

i 

If the transformation is canonical, the Kamiltonian K should fulfill a rela- 
tionship similar to (7). In other words, for the new set of variables Q and 
P we still have 

sJ{J2 P^dQi - Kdt) = . (8) 



^Here we follow the terminology of Goldstein by referring to K = K(Q,P,i), which 
is different from the Hamiltonian H = H{p,q,t) by an additive time derivative, as the 
Kamiltonian. 

^An alternative derivation has been given by G. S. S. Ludford and D. W. YanniteU, 
Am. J. Phys. 36, 231 (1968). 
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Moreover, according to the Legendre transformation, J^iPidQi — Hdt 
L{q, q, t)dt, (7) - like (8)- is equivalent to 

rt2 



6 r L{q,q,t)dt = . (9) 



In addition, (9) does not change if L is replaced by £ = L H because 

in this case 

'■^a rt2 dF{q,t), 



6 



f'cdt = 6f\L + ^^^)dt, (10) 
Jti Jti dt 



or similarly 

6 r Cdt = 6 r L{q, q, t)dt + 6F{q(^2) - 5F{q(^^^ ,ti) . (11) 

Thus, (10) and (11) differ only by constant terms whose variation is zero in 
Hamilton's principle. 

It follows that the Hamiltonian H and the Kamiltonian K are related by 
the equation ^ 

dF 

Piqi -H = P,Qi-K+ — . (12) 

The function F is the so-called generating function. It can be expressed 
as a function of any arbitrary set of independent variables. However, some 
very convenient results are obtained if F is expressed as a function of the n 
old variables and the n new ones, plus the time. The results are especially 
convenient if the n old variables are exactly the n qi - or the n pi-, and if 
the new variables are all of them the n Qi - ox the n Pi . 

Using these coordinates, the possible combinations of n old variables and n 
new variables -including t- in the generating function are^ 

Fi = Fi(Q,(?,i) (13) 

F2 = F2{P,q,t) 

Fs = Fs{Q,p,t) 

F4 = F^{P,p,t) . 



^Some authors add to the right hand side of this equation a constant multiphcative 
factor A that does not change (9). Here we use A = 1, that is, we decided to work 
with the so-caUed reduced canonical transformations, since this simpler case is sufficient 
to illustrate the structure of the canonical transformations. 

''We shall use the convention of Goldstein to denote each of the different combinations 
of the new and old variables in the generating function. 
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On the other hand, if we multiply (12) by dt we get: 



Pidqi - Hdt = PdQi - Kdt + dF . (14) 

Making the change F ^ Fi above, and recalling that dQi, dqi, and dt are 
independent variables, we get: 



Pi 



Pi 



K = H + 



dFi 

~dQi 
dFi 

dqi 

dFi 



dt 



Using now some algebraic manipulation it is possible to obtain analogous 
expressions to the previous one for the rest of the generating functions. The 
results are the following: 

F2: Q.= It m= H K = H + ^ 

F,: Q.= It q^ = -^ K = H + ^. 

In practice, one usually applies a useful theorem (see below) that allows, 
together with the definitions we gave in the introduction for canonical trans- 
formations, to solve any mechanical problem of interest Q 
Theorem 5.1 We consider a system acted by a given external force. We 
also suppose that the dynamical state of the system is determined by a set 
of variables q,p = qi,q2, ■■■,qn,Pi,P2, ■■■,Pn and that the Hamiltonian of the 
system is H = H{q,p,t). The time evolution of the variables q and p is 
given by Hamilton's eqs. 

. _ dH{q,p,t) 

Qi — 



Pi 



dpi 
dHiq,p,t) 



dqi 

If we now perform a transformation to the new variables 

Q = Q{q,p, t) ; P = P{q,p,t) 

^For an example, see the final section of this chapter. 
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and if the transformation is canonical, i.e., there exists a function F{q,p,t) 
such that for a fixed arbitrary time t = Iq we have 

dF{q,p, to) = Vidxi - ^ YidXi , 

i i 

where Xi,yi = qi,pi or pi, -Qi y Xj, Yi = Qi, Pi, or Pi, -Qi, then the equa- 
tions of motion in terms of the variables Q and P are 

. ^ dK{Q,P,t) 

dp 

dK{Q,P,t) 



dQ 



where 



dt ^ ' dt 

% 

Moreover, if the determinant of the matrix is different of zero, then 

the latter equation takes the form 

dF{x,X,t) 



K = H + 



dt 



5.2 Necessary and sufficient conditions for a trans- 
formation to be canonical 

We have already mentioned that by a canonical transformation we mean a 
transformation, which, independently of the form of the Hamiltonian, keeps 
unchanged the form of Hamilton's equations. However, one should be very 
careful with this issue because some transformations fulfill this requirement 
only for a particular Hamiltonian ^ Some authors call such transformations 
canonical transformations w.r.t. 

To illustrate this point we use the following example, given in the paper 
of J. Hurley: Let us consider a particular physical system whose Hamilto- 
nian is 

®See, for example, J. Hurley Am. J. Phys. 40, 533 (1972). 

''See, for example, R. A. Matzner and L. C. Shepley, Classical Mechanics (Prentice 
Hall, 1991). 
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and the following transformations 



It is easy to show that the Kamiltonian K given by 

2p3/2 

K = 

3m 

leads to 

P = 2pp = = -^ 
oQ 

and 

. . p P^l^ dK 

Q = q = — = = • (17) 

m m or 
On the other hand, if we choose the Hamiltonian 

2m 

then it is possible to find a Kamiltonian K for which the usage of the trans- 
formation equations (16) maintains unchanged the form of Hamilton's equa- 
tions. Thus, the equations (16) keeps unchanged the form of Hamilton's 
equations only for a particular Hamiltonian. 

It can be shown that the necessary and sufficient conditions for the canonic- 
ity of transformations of the form (3) and (4), that is, to keep unchanged the 
form of Hamilton's equations whatever the Hamiltonian, are the following 

[QuPj] = a (18) 

[P^,Pj]=0 (19) 

[Qi,Qj]=0, (20) 

where a is an arbitrary constant related to scale changes. Finally, we would 
like to make a few important comments before closing this section. First, we 
should keep in mind that Q and P are not variables defining the configuration 
of the system, i.e., they are not in general a set of generalized coordinates 
Pl To distinguish Q and P from the generalized coordinates q and p, one 



Except for the trivial case in which the canonical transformation is Q = g and P — p. 
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calls them canonical variables. In addition, the equations of motion -similar 
in form to the Hamiltonian ones- for the generalized coordinates q and p- 
that one gets for Q and P are called canonical Hamilton equations. Second, 
although we did not check here, if the transformation Q = Q(q,p,t) and 
P = P{q,p, t) is canonical, then its inverse q = q{Q, P, t) and p = p{Q, P, t) 
is also canonical 0. 



5.3 Example of application of TC 

As we already mentioned in the introduction, the main idea in performing 
a canonical transformation is to find a phase space coordinate system for 
which the form of the Hamilton eqs is maintaind whatever the Hamiltonian 
and to choose the one that makes easy the solution of the problem. We 
illustrate this important fact with the following example. 
EXAMPLE : 

The Hamiltonian of a physical system is given hy H = io'^p{q + 1)'^, where to 

is a constant. Determine g as a function of time. 

Solution: 

1. Solving the Hamilton equations for the variables q and p. Applying (1) 
and (2) to the given Hamiltonian we get 

uj'^{q + tf = q, 2uj^p{q + t) = -p . 

This system is not easy to solve. However, we can get the solution by means 
of an appropriate canonical transformation as we show in the following. 

2. Using Q = q + t, P = p. According to the theorem given in section (5.1), 
since 

'-r = » 

op 
dP 



d{-q) 

then the Kamiltonian K of the system is given by 



dF(q,p,t) do dP 



^For a proof see, for example, E. A. Desloge, Classical Mechanics, Volume 2 (John 
Wiley & Sons, 1982). 
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The form of the function F(q,p, t) can be obtained from the canonical trans- 
formation given in section 5.1 (the case corresponds to a time-dependent 
canonical transformation). Therefore, we substitute Q = q + t, P = p in 

dF{q,p, t) = pdq — PdQ , 

to get without difficulty 

F{q,p,t) = c, c= constant . 

On the other hand, 

dP 

m = ' 

dQ 
dt 

Finally, substituting these results in (21) (and also Q = q + t, P = p in H) 
we get 

K = P{lo^Q^ + 1) . 



Moreover, from (5) we find 



Q = cj'^Q^ + 1 . 



This differential equation is now easy to solve, leading to 

q = —ta,n{ujt + cj)) — t , 

CO 

where (p is an arbitrary phase. 
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6. POISSON BRACKETS 



Forward: The Poisson brackets are very useful analytical tools for the 
study of any dynamical system. Here, we define them, give some of their 
properties, and finally present several applications. 



CONTENTS: 

1. Definition and properties 

2. Poisson formulation of the equations of motion 

3. The constants of motion in Poisson formulation 
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1. Definition and properties of Poisson brackets 



If u and V are any two quantities that depend on the dynamical state of a 
system, i.e., on p and q) and possibly on time, the Poisson bracket of u and 
V w.r.t. a set of canonical variables q and p is defined as follows 

^ f du{q,p,t) dv{q,p,t) du{q,p,t) dv{q,p,t)\ 
t^'^^"?(v^^^ d^. d^. d^J- 

The Poisson brackets have the following properties (for u, v, and w arbitrary 
functions of q, p, and t; a is an arbitrary constant, and r is any of qi, pi or 

t)B 

1. [u,v] = -[v,u] 

2. [m, n] = 

3. [m, V + w\ = [u, v\ + [u, w] 
4- [u,vw\ = u[n, + [n, ujw 
5. o[u, v\ = [au, v] = [u, av] 

7. The Jacobi identity, [n, [1^,^;]] + [v, [w,u]] + [w, [u,v]] = . 

Another very important property of PBs is the content of the following 
theorem 

Theorem 6.1 If the transformation Q = Q{q,p,t), P = P{q,p,t) is a 
canonical transformation, the PB of u and v w.r.t. the variables q, p is 
equal to the PB of u and v w.r.t. the set of variables Q, P, i.e., 

/ du{q,p,t) dv{q,p,t) _ du{q,p,t) dv{q,p,t) 
^ \ dqi dpi dpi dqi 



E 



/ du{q,p,t) dv{q,p,t) _ du{q,p,t) dv{q,p,t) 
\ dQi dP^ dPi dQ^ 



As in the previous chapter, by q and p we mean q = qi, 52, Qn y P = Pi,P2, ■■■,Pn- 
^The proof of these properties can be obtained by using the definition of the PBs in 
order to express each term of these identities by means of partial derivatives of u, v, and 
w, and noticing by inspection that the resuhing equations do hold. 
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2. Poisson formulation of the equations of motion 



In the following, we outline as theorem-like statements the most important 
results on the PB formulation of the eqs of motion of the dynamical systems 



Theorem 6.2 Consider a system whose dynamical state is defined by the 
canonical variables q, p and whose dynamical behaviour is defined by the 
Hamiltonian H = H(q,p,t). Let F be an arbitrary quantity depending on 
the dynamical state of the system, i.e., on q, p, and possibly on t. The rate 
of change in time of F is given by 



where [F, H] is the PB of F and H. 

Theorem 6.3 (Poisson formulation of the eqs. of motion). Con- 
sider a system described in terms of the canonical variables q, p, and whose 
Hamiltonian is H = H{q,p, t). The motion of the system is governed in this 
case by the equations 



3. Constants of motion in Poisson's formulation 

We shall use again a theorem-like sketch of the basic results on the constants 
of motion in Poisson's formulation. These results are the following. 
Theorem 6.4 If one dynamical quantity F is not an explicit function of time 
and if the PB of F and H is zero, [F, H] = 0, then F is a constant/integral 
of motion as one can see from the theorem 6.2. 

Corollary 6.4. If the Hamiltonian is not an explicit function of time, then 
it is a constant of motion. 



^^The proofs have been omitted as being well known. See, for example, E. A. Desloge, 
Classical Mechanics, Volume 2 (John Wiley & Sons, 1982). 




qi = [qi,H] 
Pi = \PI,H] . 
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7. HAMILTON- JACOBI EQUATIONS 



Forward: It is known from the previous chapters that in principle it is 
possible to reduce the complexity of many dynamical problems by choosing 
an appropriate canonical transformation. In particular, we can try to look 
for those canonical transformations for which the Kamiltonian K is zero, a 
situation leading to the Hamilton-Jacobi equations. 

CONTENTS: 

7.1 Introduction 

7.2 Time-dependent Hamilton-Jacobi equations 

7.3 Time-independent Hamilton-Jacobi equations 

7.4 Generalization of the Hamilton-Jacobi equations 

7.5 Example of application of the Hamilton-Jacobi equations 
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7.1 Introduction 



In order to reach the goals of this chapter we need to make use of the 
following result allowing us to find the set of canonical variables for which 
the Kamiltonian takes a particular form. 

Theorem 7.1 Consider a system whose dynamical state is defined by p, q 
and whose behaviour under the action of a given force is governed by the 
Hamiltonian H = H{q,p,t). Let K = K{Q,P,t) be a known function of 
the canonical variables Q, P, and time. Then, any function F(q,Q,t) that 
satisfies the partial differential equation 



K 



dF{q,Q,t) 



dQ 

and also the condition 



H 



dF{q,Q,t) 



d^F{q,Q,t) 



dqjdQj 



dq 



+ 



dF{q,Q,t) 
dt 



is a generating function for a canonical transformation of q, p to Q, P, and 
the corresponding Kamiltonian is K = K{Q, P, t). 

In the following sections we shall use this theorem to find those canonical 
transformations whose Kamiltonian is zerc0, that leads us to the Hamilton- 
Jacobi equations. 



7.2 Time-dependent HJ equations. 



As a consequence of Theorem 7.1 and of requiring a zero Kamiltonian we 
get the following theorem 

Theorem 7.2 Consider a system of / degrees of freedom defined by the set 
of variables q, p and of Hamiltonian H = H{q,p,t). If we build the partial 
differential equation 



H 



dS{q,t) 
dq 



+ 



dSiq,t) 
dt 







(1) 



and if we are able to find a solution of the form 



5 = S{q, a, t) 



^More exactly, K [Q, -^^l^^.t] = 0. 
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where a = ai, a2, a/ is a set of constants and if in addition the solution 
satisfies the condition 

d'^S{q,a,t) 



dqida 

then q(t) can be obtained from the equations 

dS{q,a,t) 



dai 



(2) 



where /? = /3i, /32, /3/ is a set of constants. The set of equations (2) provide 

us with / algebraic equations in the / unknown variables qi,q2, ■■■,qf- The 
values of the constants a and /? are determined by the boundary conditions. 
Moreover, if q(t) is given it is possible to find p(t) starting from 



Pi 



dS{q, a, t) 
dqi 



(3) 



The partial differential equation (1) is called the time- dependent Hamilton- 
Jacobi equation. The function S{q, a, t) is known as Hamilton's principal 
function. 

To achieve a better meaning of the theorem, as well as of the constants a 

and P, we proceed with its proof. 

Proof of the Theorem 7.2. According to the Theorem 7.1, any function 
F{q,Q,t) satisfying the partial differential equation 



H 



dF{q,Q,t) , 
q, „ : E 



and also the condition 



dq 

d^F{q,Q,t) 



I dF{q,Q,t) _^ 



dt 



^0 



dqjdQj 

should be a generating function of a set of canonical variables Q, P for which 
the Kamiltonian K is zero, i.e., K[Q,P,t) = 0. The function 

F{q, Q, t) = [S{q, a, t)U=Q = S{q, Q, t) 

belongs to this class. Then S{q, Q,t) is the generating function for a canon- 
ical transformation leading to the new set of canonical variables Q, P, for 
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which the Kamiltonian K is identically zero. The transformation equations 
associated to S(q,Q,t) are 



Pi 
Pi 



dS{q,Q,t) 

dS{q,Q,t) 
dQi 



and because K{Q,P,t) = 0, the equations of motion are 

A _ dK(Q,P,t) _ 

— p)p. — u 

Pi 

Prom these equations we infer that 



_ dK(Q,P,t) _ 



Qi = ai 
Pi = -P, 



(4) 
(5) 



(6) 
(7) 



where ai and Pi are constants. The choice of the negative sign for j3 in (7) 
is only a convention. If now we substitute the equations (6) and (7) in (5) 
we get 



-Pi = - 



dS{q,Q,t) 
dQi 



Q=a 



dS{q, a, t) 
dai 



which reduces to (2). If, in addition, we substitute (6) in (4) we get (2). 
This complets the proof. 



7.3 Time-independent HJ equations 



If the Hamiltonian does not depend explicitly on time, we can partially solve 
the time-dependent Hamilton-Jacobi equation. This result can be spelled 
out as the following theorem 

Theorem 7.3 Consider a system of / degrees of freedom defined in terms 
of q, p, and whose behaviour under a given force is governed by the time- 
independent Hamiltonian H{q,p). 
If we build the partial differential equation 



H 



dWjqY 
dq 



E, 



(8) 
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where £J is a constant whose value for a particular set of conditions is equal to 
the value of the integral of motion H{q,p) for the given boundary conditions, 
and if we can find a solution to this equation of the form 

W = W{q,a) , 

where a = ai,a2) •••,«/ is a set of constants that explicitly or implicitly 
include the constant E, i.e., E = E{a), and if the solution satisfies the 
condition 

d'^W{q,a) 



dqidaj 

then the equations of motion are given by 

dS{q, a, t) 



^0, 



dai 

where 



= A (9) 



S{q, a, t) = W{q, a) - E{a)t 

and (5 = (3i, P2, Pf is a set of constants. The set of equations (9) pro- 
vide / algebraic equations in the / unknown variables gi, 52, The 
values of the constants a and /3 are determined by the boundary conditions. 
The partial differential equation (8) is the time-independent Hamilton- Jacobi 
equation, and the function W{q,a) is known as the characteristic Hamilton 
function. 



7.4 Generalization of the HJ equations 

The Hamilton-Jacobi equation can be generalized according to the follow- 
ing theorem allowing sometimes the simplification of some Hamilton-Jacobi 
problems. 

Theorem 7.4 Consider a system of / degrees of freedom whose dynamics 
is defined by x, y, where Xi, yi = qi,Pi or pi, —qi, and whose behaviour under 
the action of a given force is governed by the Hamiltonian H{x,y,t). If we 
write the partial differential equation 



H 



' dS{x,t) , 
ox 



+ 



dS{x,t) 
dt 







and if we can find a solution of this equation of the form 



S = S{x,a,t) 
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where a = ai, 02, a/ is a set of constants and in addition the solution 
satisfies the condition 

d'^S{x,a,t) 



dxjdaj 



then the laws of motion of the system can be obtained from the equations 

dS{x, a, t) 



dxi 
dS{x, a, t) 
dai 



(10) 

(11) 



where /? = /?2, /3/ is a set of constants. 



7.5 Example of application of the HJ equations 

We shall solve the problem of the one-dimensional harmonic oscillator of 
mass m, using the Hamilton-Jacobi method. 
We know that the Hamiltonian of the system is 

H = ^ + ^ (12) 
According to Theorem 7.2 the Hamilton-Jacobi equation for the system is 



2^ Vdq) ^~^~dt ~ 



(13) 



2m\dq ) 2 dt 
Making each side of the previous equation equal to a, we find 

2 ;_2 



Wc assume a solution of (13) of the form F = Fi(q) + F2{t). Therefore, (13) 
converts to 

1 fdFi\\ kq^ dF2 ^^^^ 



(15) 
(16) 

(17) 
(18) 



J_ /dF\ 
2m \ dq 



kg' 

+ - = a 



dF2 
dt 



= —a 



For zero constants of integration, the solutions are 



Fi = J ^2m{a-^)dq 
F2 = 



-at 
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Thus, the generating function F is 




According to (2), q{t) is obtained starting from 




and effecting the integral we get 

sm-^{q^k/2a) = t + f3 , (22) 
from which q is finally obtained in the form 

q = sin k/m{t + P) . (23) 

In addition, we can give a physical interpretation to the constant a according 
to the following argument. 

The factor should correspond to the amplitude A of the oscillator. On 
the other hand, the total energy E oi a, one-dimensional harmonic oscillator 
of amplitude A is given by 




In other words, a is physically the total energy E of the one-dimensional 
harmonic oscillator. 
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8. ACTION-ANGLE VARIABLES 



Forward: The Hamilton-Jacobi equation provides a link to going from a 
set of canonical variables q, pto a second set Q, P, where both are constants 
of motion. 

In this chapter, we briefly present another important procedure by which 
one goes from an initial pair of canonical variables to a final one, where not 
both variables are simultaneously constants of motion. 

CONTENS: 

8.1 Separable systems 

8.2 Cyclic systems 

8.3 Action-angle variables 

8.4 Motion in action-angle variables 

8.5 Importance of action-angle variables 

8.6 Example: the harmonic oscillator 
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8.1 Separable systems 



Separable systems are those ones for which the Hamiltonian is not an ex- 
pHcit function of time, i.e., 

H = H{q,p) , 

allowing, in addition, to find a solution of the time-independent Hamilton- 
Jacobi of the form 

W{q,a) = Y,Wi{qi,a) . 

i 

8.2 Cyclic systems 

We know that the dynamical state of a system is characterized by a set of 
generalized coordinates q = qi,q2, ---jqf and momenta p = pi,p2, ■.■,pf. A 
system in motion will describe an orbit in the phase space q, p At the same 
time, there is an orbit in each of the subspaces qi,Pi. In every plane qi,Pi, 
the orbit can be represented by an equation of the form pi = Pi{qi) or a pair 
of equations pi = Pi{t), qi = qi{t). If for each value of i, the orbit Pi = Pi{qi) 
is a closed curve in the plane qi — Pi, then we say that the system is cyclic. 
In the figure 8 we show the two possibilities for a system to be cyclic. In 
8.1a, the system is cyclic because q^ oscillates between the limits defined by 
qi = a and q^ = b, whereas in figure 8.16, the system is cyclic because qi 
moves from qi = a to qi = b, and repeats the same motion afterwards. 




b 



Pig. 8 .1 

At this point, it is worthwhile to make two helpful remarks. 

Remark 1: The cyclic term has been introduced only for simplifying 
the notation in the next sections. One should not interpret this term as if 
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the system is cyclic in each subspace qi,pi. The system should come back 
to its initial state only in the global space q,p. 

Remark 2: If the cyclic system has only one degree of freedom, the time 
required by the system to accomplish the cycle in q — p is constant; therefore 
the motion in the space q—p will be periodic in time. If the system has 
more degrees of freedom, then, in general, the time required for a particular 
cycle in one of the subspaces qi,Pi will not be a constant, but will depend on 
the motion of the other coordinates. As a result, the motion in the subspace 
qi,Pi will not be periodic in time. One should be careful with this point, 
since not all the motions in the subspaces qi,Pi are periodic. 

8.3 Action-angle variables 

We consider now a cyclic system of / degrees of freedom, whose dynamical 
state is characterized by the canonical set q, p. Let H{q,p) be the Hamilto- 
nian of the system and let 



W{q,a)=J2Wiiqi 



a) 



(where a = ai,a2, ---jaf are constants) be a solution of the time- independent 
Hamilton-Jacobi equation 

dW 

Let J = Ji, J2, J/ be the set of constants defined by the equations 

JiW = f Q^, dqi , (1) 

where the integral is along a complete cycle in the variable qi. If we use the 
function 

W{q,a) = W[q,aiJ)] 

= 5^W-4g,,a(J)] 

i 
i 

as the generating function of a canonical transformation of q, p to a new set 
of coordinates w and momenta J, i.e., if we define the variables w and J by 
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the transformation equations 



^ dW{q,a) ^ dWi{qi,J) 
dqi dqi 

dWiq,J) 

= -^j— , (3) 

then the new coordinates ■wi,W2, ■■■,Wf are called angle variables, and the 
new momenta Ji, J2, J/ are called action variables. 
From (2) we get 

^W^[qi,J{a)] dWi{qi,a) 

Pi{qi,a) = = . (4) 

dqi dqi 

Substituting (4) in (1) one gets 

Ji{a) = f pi{qi,a)dqi . (5) 



The equation pi = Pi{qi,a) gives the projected orbit p = p{q) on the subspace 
Pi,qi- The integral in the right hand side of the equation (5) is thus the area 
bordered by the closed orbit, or beneath the orbit, as shown in figure 8.1. 
Thus, the function Jj(a) has a geometric interpretation as the area covered 
in the subspace qi,Pi during a complete cycle in the subspace. This area 
depends on the constants a or equivalently on the initial conditions and can 
be arbitrary |^. 



8.4 Motion in terms of action-angle variables 

We give the following theorem-like statement for the motion of a system in 
terms of action-angle variables. 
Theorem 8.4 

Consider a separable cyclic system of / degrees of freedom whose motion 
is described by the variables q,p = qi,q2, ■■■,qf,Pi,P2, together with 

the Hamiltonian H{q,p). If we transform the motion to the action-angle 
variables J,w, then the Hamiltonian H is a function of J alone, i.e., 

H = H{J) 

Historically, the first intents to pass from the classical mechanics to quantum mechan- 
ics was related to the assumption that the value of J; could be only a multiple of h/2-n, 
where h is Planck's constant. 
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and the equations of motion will be 



Ji = 

Wi = 



vd + 1 



where 7^ y arc constants determined by the initial conditions, whereas 
the Vi are constants known as the. frequencies of the system being defined as 
follows 



-dHjjy 

_ dJi _ 



8.5 Importance of the action-angle variables 

The importance of the action-angle variables resides in providing a powerful 
technique for directly getting the frequencies of periodic motions without 
solving for the equations of motion of the system. 

This important conclusion can be derived through the following argument. 
Consider the change of w when q describes a complete cycle 

/ 9w , 
Aw = (p -T—dq 
J oq 

On the other hand, we know that 

dW 



and therefore 




= 1 . 



This result shows that w changes by a unity when q varies within a complete 
period. 

Prom the relationship 

W = L't + (f), 
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we infer that in a period r 



Alt; = 1 

= VT 



This means that we can identify the constant u with the inverse of the period 

1 

V = — . 

T 

8.6 Example: The simple harmonic oscillator 

Using the action-angle formahsm prove that the frequency v of the simple 
one-dimensional harmonic oscillator is given hy v = ^fkJrajl'K. 
Since is a constant of motion, the orbit in the space q — p'x^ given 

^ + ^ = ^, 
2m 2 

where E is the energy. This is the equation of an ellipse of semiaxes 

and y^2E/k. The area enclosed by the ellipse is equal to the action J. 

Therefore, 

/TT—^ l^E RE^ 
J = nV2rnE\ — = 2ttJ—E . 

V A; V A; 

It follows that 



H{J) = E = ^J. 

27r 



Thus, the frequency will be 



_ dH{J) _ ^Jkjm 
^~ dJ ~ 27r 
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9. CANONICAL PERTURBATION THEORY 



Forward: The great majority of problems that we have to solve in Physics 
are not exactly solvable. Because of this and taking into account that we 
live in the epoch of computers the last decades have seen a lot of progress in 
developing techniques leading to approximate solutions. The perturbation 
method is used for not exactly solvable Hamiltonian problems when the 
Hamiltonian differs slightly from an exactly solvable one. The difference 
between the two Hamiltonians is known as the perturbation Hamiltonian. 
All perturbation methods are based on the smallness of the latter with 
respect to both Hamiltonians. 



CONTENTS: 

9.1 Time-dependent perturbation theory (with two examples) 

9.2 Time-independent perturbation theory (with an example) 
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9.1 Time-dependent perturbation theory 



The most appropriate formulation of classical mechanics for the develop- 
ment of perturbation methods is the Hamilton-Jacobi approach. Thus, we 
take Ho{p,q,t) as the Hamiltonian corresponding to tlic solvable (unper- 
turbed) problem and consider the Hamilton principal function S{q, ao,t), as 
the generating function of a canonical transformation of {p, q) to the new 
canonical pair (ao,/3o) for which the new Hamiltonian (or Kamiltonian) Kq 
of the unperturbed system is zero 



8^ r)"? 

-+m-.,,t) = K„=o. (1) 

This is the Hamilton-Jacobi equation, where we used p = dS/dq. The whole 
set of new canonical coordinates (q:o,/3o) are constant in the unperturbed 
case because Kq = and: 



^ f^. (2) 
We now consider the Hamiltonian of the perturbed system written as follows: 



H{q,p,t)=Ho{q,p,t)+AH{q,p,ty, {AH ^ Ho). (3) 

Although (ao,/3o) a-rc still canonical coordinates (since the transformation 
generated by S is independent of the particular form of the Hamiltonian), 
they will not be constant and the Kamiltonian K of the perturbed system 
will not be zero. In order not to forget that in the perturbed system the 
transformed coordinates arc not constant, we denote them by a and (3 in- 
stead of cto and /3o, which are the corresponding constants in the unperturbed 
system. Thus, for the perturbed system we have: 

(99 dS 

K{a, p,t) = H + — = iHo + + AH = AH{a, /3, t). (4) 

The equations of motion for the transformed variables in the perturbed 
system will be: 
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Oli = — 



Pi = 



dA H{a,l3,t) 

dAH{a,(3,t) 
dai 



(5) 



where i = 1, 2, n and n is the number of degrees of freedom of the system. 
These are rigorous equations. If the system of 2n equations could be solved 
for ai and (3i as functions of time, the transformation equation (p, g) — > (a, (5) 
would give pi and qi as functions of time and the problem would be solved. 
However, the exact solution of the equations (5) is not easier to get w.r.t. 
the original equations. From (5) we see that when a and [3 are not constant 
their variation in time is slow if we assume that AiJ changes infinitesimally 
w.r.t. a y p. A first approximation for the temporal variations of (a, /3) 
can be obtained by substituting a y /3 in the second terms of (5) by their 
constant unperturbed values 



an 



dAH{a,f3,t) 

m 

dAH{a,p,t) 
dai 



(6) 



where an and Pn are the first-order solutions, i.e., in the first power of the 
perturbation to a^ and Pi, and the vertical bars with zero subindices show 
that after the derivation one should substitute a and P by their constant 
unperturbed values. Once this is done, the equations (6) can be integrated 
leading to ai and Pi as functions of time (in the first order). Next, using 
the equations of transformation one can get p and q as functions of time in 
the same first-order approximation. The second-order approximation can 
be obtained now by substituting in the second terms of (6) the first-order 
approximation of a and P w.r.t. time. In general, the perturbation solution 
of order N is obtained by integrating the following equations 



OliN 
PiN 



dAH{a,P,t) 

dPi 
dAH{a,P,t) 
dai 



N-l 



N-1 



(7) 



Example 1 



Let us consider the simple case of a free particle that next is subject to 
a harmonic perturbation. Although this example is trivial can be used to 
illustrate the aforementioned procedure. The unperturbed Hamiltonian is 
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(8) 



Since Hq ^ Ho{x), x is a cyclic variable and p = ao is a constant of motion 
in the unperturbed system. Recalling that p = dS/dx, we substitute in (1): 



1 fdSy dS 



Since the system is conservative, it is convenient to consider the principal 
function of the form 



S = Six)+Fit). (10) 



This type of separation of variables is quite useful when the Hamiltonian 
does not depend explicitly on time. Then, one writes F{t) = —Et, where E 
is the total energy of the system |^. Putting (10) in (9) we obtain 



2m \dx J 



] = E, S = V2mEx = aox. 



(11) 



Substituting (11) in (10), together with the fact that in this case the Hamil- 
tonian is equal to the energy, we can write the principal function of Hamilton 
as follows 



5 = aox-—. (12) 
2m 



If the transformed momentum is qq, the transformed coordinate (which is 
also constant in the unperturbed system) is 



dS aot 

Po = T, — = x 

oao m 



Therefore, the transformation generated by S is given by the following 

eqs 



'See, M.R. Spiegel, Theoretical Mechanics, pp. 315, 316. 
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P = "0, 




(13) 



m 



They represent the solution for the motion of free particle. What we have 
given up to now is only the procedure to obtain the equations of motion using 
the Hamilton-Jacobi formulation. Now, we introduce the perturbation 



or, in terms of the transformed coordinates, using (13) 



Notice that we have renounced at the subindices for the transformed co- 
ordinates since we already study the perturbed system. 

Substituting (15) in (5) we get 



As one may expect, these equations have an exact solution of the harmonic 
type. To be sure of that we perform the time derivative of the first equa- 
tion allowing us to conclude that a has a simple harmonic variation. The 
same holds for x as a consequence of the transformations (13), which are 
invariant in form in the perturbed system (up to the subindices of the trans- 
formed coordinates). However, we are interested to illustrate the pertur- 
bation method, so that we consider that k (the elastic constant) is a small 
parameter. We seek approximate solutions in different perturbative orders, 
without missing the fact that the transformed variables (a,/?) in the per- 
turbed system are not constants of motion. In other words, even though 
(a, P) contain information on the unperturbed system, the effect of the per- 
turbation is to make these parameters varying in time. 

The first-order perturbation is obtained in general as given by (6). Thus, we 




(15) 




(16) 
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have to substitute a and (3 by their unperturbed values in the second terms 
of (16). To simphfy, we take x{t = 0) = and therefore /3o = 0, leading to 



ai = —uj^aot, 

k = ao^, (17) 
m 



which integrated reads 



The first-order solutions for x and p are obtained by putting ai and /?i in 
the transformation eqs. (13), where from 



X = — , 

mu \ b I 

P = «ofl-^^V (19) 



To obtain the approximate solution in the second perturbative order we have 

to find d2 and (32, as was indicated in (7), by substituting in the second terms 
of (16), Q and /3 by ai and /3i as given in (18). Integrating d2 and (32 and 
using again the transformation eqs. (13), we get the second-order solutions 
for X and p: 



P = «o 1-^ + ^ . (20) 



In the limit in which the perturbation order tends to infinity, we obtain 
the expected solutions compatible with the initial conditions 




X — ^ sina;t, p—>aocosu!t. (21) 



The transformed variables (a, (3) contain information about the unperturbed 
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orbit parameters. For example, if we consider as unperturbed system that 
corresponding to tlie Kepler problem, a convenient coordinate pair (a, (3) 
could be the (J, S) variables, which are the action and the phase angle of the 
angle w, respectively (remember that w = vt + S, where v is the frequency). 
These variables are related to the set of orbital parameters such as semimajor 
axis, eccentricity, inclination, and so on. 

The effect of the perturbation is to produce a time variation of all these 
parameters. If the perturbation is small, the variation of the parameters 
during a period of the unperturbed motion will also be small. In such a 
case, for small time intervals, the system moves along the so-called osculat- 
ing orbit, having the same functional form as the orbit of the unperturbed 
system; the difference is that the parameters of the osculating curve vary in 
time. 

The osculating parameters can vary in two ways 

• Periodic variation: if the parameter comes back to its initial value 
after a time interval that in first approximation is usually the unper- 
turbed period. These periodic effects of the perturbation do not alter 
the mean values of the parameters. Therefore, the trajectory is quite 
similar to the unperturbed one. These effects can be eliminated by 
taking the average of the perturbations in a period of the unperturbed 
motion. 

• Secular variation: At the end of every succesive orbital period there 
is a net increment of the value of the parameter. Therefore, at the 
end of many periods, the orbital parameters can be very different of 
their unperturbed values. The instantaneous value of the variation of 
a parameter, for example the frequency, is seldom of interest, because 
its variation is very small in almost all cases in which the perturbation 
formalism works. (This variation is so small that it is practically im- 
possible to detect it in a single orbital period. This is why the secular 
variation is measured after at least several periods.) 



Example 2 

From the theory of the Kepler two-body problem it is known that if we add 
a potential the orbit of the motion of negative energy is a rotating 

ellipse whose periapsis is precessing. In this example, we find the precession 
velocity for a more general perturbation 
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where n > 2 is an entire number, and h is such that the second term of 
the potential is a small perturbation of the first one. The perturbative 
Hamiltonian is 



^H = -^. (23) 



In the unperturbed problem, the angular position of the periapsis in the 
orbit plane is given by the constant uj = 2'kw2- In the perturbed case, we 
have 



. ^ dAH dAH 

^ = ^^^ = ^^ ^''^ 

where we have used J2 = ^ttI. Moreover, J2 and W2 are two of the five 
integrals of motion that can be obtained when one uses the action-angle 
variables to solve the Kepler problem. 

We need to know the mean of cj in a period r of the unperturbed orbit 

, . , 1 r dAH , d fl r ^ \ d{AH) , , 

The temporal mean of the unperturbed Hamiltonian is 

= (26) 

On the other hand, since I = mr'^{d6/dt), we can get dt and plunge it in 
(27). This leads to 



{AH) = - — 



It Jo r" ^ 

= -17(-^j X + (27) 

where is a constant phase, e is the eccentricity, and where we expressed r 
as a function of 9 using the general equation of the orbit with the origin in 
one of the focal points of the corresponding conic 
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1 

r 



mk 



1 + e cos {6 — rj) 



(28) 



For n = 2: 



It 

2'Kmh 

i^) = (29) 



For n = 3: 



(Ail) 

(a;) 



2'Km?hk 



(30) 



The latter case, n = 3, is of particular importance because the theory of 
General Relativity predicts a correction of the Newtonian motion precisely of 
order. This prediction is related to the famous problem of the precession 
of the orbit of Mercury. Substituting the appropriate values of the period, 
mass, great semiaxis, that is included in h, and so on, (30) predicts a mean 
precession velocity of 



(uj) = 42.98 arcsec. /century . 



The mean value is by far larger than the aforementioned one (by a factor 
larger than one hundred). But before making any comparison one should 
eliminate from the mean value the contributions due to the following factors: 
a) the effect known as the precession of the equinoxes (the motion of the ref- 
erence point of longitudes w.r.t. the Milky Galaxy), b) the perturbations of 
the Mercury orbit due to the interaction with the other planets. Once these 
are eliminated, of which the first is the most significant, one may hopefully 
obtain the contribution of the relativistic effect. In 1973, this contribution 
has been estimated as 41.4ib0.9 arcsec. /century. This is consistent with the 
prediction given by (30). 

9.2 Time-independent perturbation theory 

While in the time-dependent perturbation theory one seeks the time depen- 
dence of the parameters of the unperturbed system initially considered as 
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constant, the aim of the time- independent approach is to find the constant 
quantities of the perturbed system. This theory can be apphed only to con- 
servative periodic systems (both in the perturbed and unperturbed state). 
For example, it can be applied to planetary motion when one introduces any 
type of conservative perturbation to the Kepler problem, in which case it is 
known as the von Zeipel or Poincare method. 

Here, we consider here the case of a periodic system of one degree of freedom 
and time-independent Hamiltonian of the form 



H = H{p,q,X), (31) 



where A is a small constant specifying the strength of the perturbation. We 
assume that 



Hoip,q)=H{p,q,0) (32) 

corresponds to a system that has an exact (closed-form) unperturbed solu- 
tion in the action-angle variables (Jo,wo), i.e.. 



Hoip,q) = Kq{Jq 



vq = wo = ; {wq = uot + 60). (33) 



Since the canonical transformation from {p,q) to (Jq,wq) is independent of 
the particular form of the Hamiltonian, the perturbed Hamiltonian H(p, q, A) 
can be written as H{Jo,wo,X). Due to the fact that the perturbed Hamil- 
tonian depends on wq, Jq, it is not constant any more. On the other hand, 
in principle, one can get new action-angle variables (J, w) that may be more 
appropriate for the perturbed system, such as 



H{p,q,X) = E{J,X) 
dl 

~dj 

dw 



. = - = ^ (34) 
dE 

J = — — — = ; {J = constant). 



Since the transformation connecting (p, g) to {Jq,wq) is known, we have 
to find the canonical transformation 5 connecting {Jq,ivq) to {J,w). If we 
assume that A is small the transformation we look for should not differ much 
from the identity transformation. Thus, we write the following expansion 
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S = S{wo, J, A) = So{wo, J) + XSiiwo, J) + X''S2{wo, J) + 



(35) 



For A = we ask S to provide an identity transformation that leads to 

^0 = woJ . (36) 
The canonical transformations generated by 5" read 



dS BSi. ^. ,2^52, 

w = — = wo + \—{wo,J) + X —{wo,J) + ... . 

Jo = — = J + A-^ u;o,J «^o,^ +••■ ■ 37 

OWq OWq OWq 

Due to the fact that wq is an angle variable of the unperturbed system we 
know that Awq = 1 over a cycle. On the other hand, we know that the 
canonical transformations have the property to conserve the phase space 
volume. Therefore, we can write: 



J = J> pdq = j> Jodwo . (38) 

Integrating the second equation of (37) along an orbit of the perturbed 
system, we get 

^ Jodwo = j) JdvjQ + ^ A" j) ^^dwQ, (39) 
and substituting (39) in (38) leads to 

J = JAwo + ^ A" / 

n=l 

Since Awq = 1, one gets 



^dwo. (40) 
own 



i:A"f||.»„ = 0. (41) 



n=l 



or 
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^dwo = 0. (42) 
own 



Moreover, the Hamiltonian can be expanded in A as a function of wq and 
Jo: 



H{wo, Jo, A) = Ko(Jo) + XKi{wo, Jo) + X'^K2iwo, Jq) + ... , (43) 



where the Ki are known because if is a known function of wq and Jq for a 
given A. On the other hand, one can write 



Hip,q,X) = H{wo,Jo,X) 

= E{J,X), (44) 

which is the expression for the energy in the new action-angle coordinates 
(where J will be constant and w will be a linear function of time). 

E can also be expanded in powers of A: 

E{J, A) = EoiJ) + XEi{J) + X^E2{J) + ... . (45) 



Taking into account (44) wc can obtain equalities for the coefficients of the 
same powers of A in (43) and (45). However, these expressions for the energy 
depend on two different sets of variables. To solve this issue we express Hq 
in terms of J by writing a Taylor expansion of H{wo, Jq, A) w.r.t. Jq in the 
infinitesimal neighborhood of J: 



H{wo, Jo, A) = H{wo, J, A) + (Jo - J)|^ + ^'"^ > '^^^ + ... , (46) 



The derivatives of this Taylor expansion, which in fact are the derivatives 
w.r.t. Jo calculated for Jq = J, can also be written as derivatives w.r.t. J, 
once we substitute Jq by J in Hq{ Jq). In the previous equation, all the terms 
containing Jq should be rewritten in terms of J by using the transformation 
defined by (37) connecting the coordinates {Jo,wo) and {J,w). Thus, from 
the second equation in (37) we obtain (Jq — J), which introduced in (46) 
gives: 
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Hi.., A) ^ Hi... A) . § (a^ . A^ll ....). 10A^ (^) ' . 0(A», ,47) 

Next, we can make use of (43) to write H{wo, J, X) = H{wo, Jq, X) \jq=j, 
that we put in (47) in order to get 



H{wo,Jo,X) = Ko{J) + XKi{wo,J) + X^K2{wo,J) + ... 



dwQ 



dj 



dJ 



+ X 



dKo{J)dS2 ^ Id^Kp (dS^y ^ 



dJ dwQ ' 2 dJ'^ \dwo 
E{J,X) 

Eo{J) + XEi{J) + X^E2{J) + .... 



(48) 



Now, we can solve the problem in terms of the coefficients Ei{ J) allowing us 
to calculate the frequency of the perturbed motion in various perturbation 
orders. Since the expansion of the terms of Ei does not imply a dependence 
on wq, then the occurrence of wq in (48) is artificial. The Ki{wo,J) of 
(48) are known functions, whereas the Si{wo, J) and Ei{J) are the unknown 
quantities. 

Making equal the corresponding powers of A, we get 



Eo{J) = Ko{J) 

dSi dKo{J) 



Ei{J) = Ki{wo,J) + 
E2{J) = K2iwo,J) + 



dwQ dJ 

dSi dKi{wo,J) 

dwQ dJ 



1 /dS.y d^KojJ) dS2 dKo{J) 

2 \dwoJ ^ dwo dJ ' ^ ^ 



We can sec that to obtain Ei we need to know not only Ki but also ^i. 
Moreover, we should keep in mind that the E^ are constant, being functions 
of J only, which is a constant of motion. We also have to notice that dKo / dJ 
does not depend on wq (since Kq = Kq{J) = Kq{Jq) | jo=j)- Averaging over 
wq on both sides of the second equation in (49), one gets 
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(El) 
{Ki) 



dJ ^dwo' 



(50) 



But we have already seen that (dSi/dwo) = §{dSi/dwo)dwo = 0. Therefore, 



El = {El) = {Ki) . (51) 

Introducing (51) in the left hand side of the second equation of (49) we get 
(dSi/dwo) as follows 

dS^ ^ {Ki)-Ki ^ 
dwQ vo{J) 

where we used uq = dK^/dJ. 

The solution for Si can now be found by a direct integration. In general, 
once we assume that we already have En-i), the procedure to obtain En is 
the following 

• Perform the average on both sides of the nth equation of (49). 

• Introduce the obtained mean value of En, in the complete equation for 
En given by (49) (the one before averaging). 

• The only remaining unknown Sn can be obtained by integrating the 
following relationship 

dS 

—-^ = known function of Wq and J . 

OWq 

• Substitute Sn in the complete equation for En- 



Once all this has been done, the procedure can be repeated for n + 1. 

As one can see, getting the energy at a particular order n is possible if and 
only if one has obtained the explicit form of Sn-i- On the other hand, Sn 
can be obtained only when En has been already found. 

The time-independent perturbation theory is very similar to the Rayleigh- 
Schroedinger perturbation scheme in wave mechanics, where one can get 
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En only if the wavefunction is known at the n — 1 order. Moreover, the 
wavefunction of order n can be found only if En has been calculated. 
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10. ADIABATIC INVARIANTS 



Forward: An adiabatic invariant is a function of the parameters and con- 
stants of motion of a system, which remains almost constant in the hmit in 
which the parameters change infinitesimally in time, even though in the end 
they may change by large amounts. 



CONTENTS: 

10.1 BRIEF HISTORY 
10.1 GENERALITIES 
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10.1 BRIEF HISTORY 



The notion of adiabatic invariance goes back to the early years of quantum 
mechanics. Around 1910 people studying the emission and absorption of 
radiation noticed that the atoms could live in quasi-stable states in which 
their energy was alsmost constant. At the 1911 Solvay Congress, the problem 
of adiabatic invariance became widely known due to Einstein that draw the 
attention of the physics community to the adiabatic invariant Ejv oi the one- 
dimensional pendulum of slowly varying length. He suggested that similar 
invariants could occur for atomic systems setting their stability limits. Later, 
Ehrenfest was able to find such adiabatic invariants and their employment 
led to the first quantum approach of Bohr and Sommerfeld. 

The method of adiabatic invariants resurged after several decades in the area 
of magnetospheric physics. The Scandinavian scientists were especially in- 
terested in boreal aurora phenomena, i.e., the motion of electrons and ions 
in the terrestrial magnetosphere. One of them, H. Alfven, showed in his 
book Cosmic Electrodynamics that under appropriate conditions a certain 
combination of dynamical parameters of the charged particles remains con- 
stant in the first order . Apparently, Alfven did not realized that he found 
an adiabatic invariant. This was pointed out by L. Landau and E. Lifshitz 
who discussed in detail this issue. 

10.2 GENERALITIES 

We first show in a simple way which is the adiabatic invariant in the case of 
the one-dimensional harmonic oscillator. The employed method will be to 
prove that in the limit of the infinitesimally slow variation of the parameters, 
the adiabatic invariant of the one-dimensional system goes to a quantity that 
is exactly conserved in a corresponding two-dimensional system. 

Consider a particle at the end of a rope of negligible mass that rotates uni- 
formly on a table. Let a be the radius of the circle, m the mass of the 
particle, and lu the angular frequency. If the origin is in the center of the 
circle the projection of the motion on the x axis corresponds to the motion 
of a simple harmonic oscillator. The motion continues to be harmonic even 
when the rope is slowly made shorter by pulling it through a small hole 
drilled at the origin. The problem is to find the invariant quantity in this 
case. It is by far more difficult to get the answer for the circular motion 
in comparison with the corresponding projected oscillator. Since there are 
only central forces even when the rope is slowly shortened the conserved 
quantity is the angular momentum I. Therefore, in the slow limit one can 
write 
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(1) 



From this equation one can see that although I is an exact invariant mcoa^ 
is only a slow limit, i.e., adiabatic invariant. More exactly, it is only in the 
slow limit that we can claim that the two quantities are equal and we are 
allowed to think of the right hand side as an invariant. In the slow limit, 
the projected onc-dimcnsional oscillator has a slowly varying amplitude, a 
frequency v = oj/Itt, and a total energy E = ^muo^a^. Therefore, the adia- 
batic invar iance of the right hand side of (1) implies the adiabatic invar iance 
of the quotient E/u for the projected oscillator. 

A short analytic proof of this fact is the following. One writes the two- 
dimensional Hamiltonian 

H={2m)-\pl+pl) + ^muj\x' + y^) . (2) 

Hamilton's equations lead to 

X + oJ^x = o, y + oJ^y = . (3) 

Changing to polar coordinates we note that / = mOr'^. The previous Hamil- 
ton equations can be written as a single complex equation z + uP'z = 0, 
where z = x + iy. We consider a solution of the form z = aexpi{Lot + b), 
where a, b, and lo are real functions that arc slowly varying in time. The real 
and imaginary parts are x = acos{u)t + b) and y = asm{ujt + b), respectively. 
Since co, a, and b vary slowly, / can be written in a very good approximation 
as I = mua^. The energy of the projected oscillator \s E = ^mx^ + ^mu^x^. 
In the slow varying limit of the parameters, we can substitute x in £' by the 
cosine function to get E = ^m.uj'^a!^. Thus, if muoa?' is an adiabatic invariant 
then E/u) is an adiabatic invariant for the one-dimensional oscillator. 

We now clarify the meaning of the term almost constant in the definition of 
the adiabatic invariant in the forward to this chapter. 

Consider a system of one degree of freedom, which initially is consevative 
and periodic and depends on an initially constant parameter a. The slow 
variation of this parameter due, for example, to a low amplitude perturba- 
tion does not alter the periodic nature of the motion. By a slow variation 
we mean the one for which a varies slowly during a period r of the motion: 

T{da/dt) < a . (4) 

However, even when the variations of a are small during a given period, after 
a sufficently long time the motion may display large changes. 
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When the parameter a is constant, the system wih be described by the 
action-angle variables (u^o, Jo) and the Hamiltonian H = H{Jq, a). Assume 
now that the generating function of the transformation {q,p) — {wq, Jq) is 
of the form W*{q, wo,a). 

When a is allowed to vary in time, (u^o, Jo) will still be valid canonical 
variables, but W* turns into a function of time through a. Then, neither 
Jo will be a constant nor wq a linear function of time. The appropriate 
Hamiltonian will be 



dW 

K{wo,Jo,a) = H{Jo,a) + 



dt 

dW* 

= H{Jo,a) + a-^ . (5) 

The second term in (5) can be seen as a perturbation. Then, the temporal 
dependence of Jq is given by 

Jo = = -d— (^] (6) 

dwo dwQ \ da J ' 

Proceeding similarly to the time-dependent perturbation theory, we seek 
the variation to first-order of the mean value of Jo during the period of the 
unperturbed motion. Since a varies slowly, we can think of it as constant 
during this interval. Thus, we can write 

1 /■ d /dW*\ , d f d /dW*\ , 2 N 

One can prove that W* is a periodic function of wo, and therefore, can be 
written, together with its derivative w.r.t. a, as a Fourier series 

^ = EA,e'^''^-'^. (8) 
Substituting (8) in (7) we get 

(Jo) = -- / ^27rifcAfce2"^'="'°di + 0(d2,d). (9) 
■^'^ k 

Since the integrant does not contain any constant term the integral is zero. 
This leads to 
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(Jo) =0 + O(d2,a). 



(10) 



Thus, (Jo) will not display secular variations in the first order, i.e., in d, 
which is one of the basic properties of the adiabatic invariance. In this way, 
the term almost constant in the definition of the adiabatic invariant should 
be interpreted as constant in the first order. 

Further reading 

• L. Parker, Adiabatic invariance in simple harmonic motion, Am. J. 
Phys. 39 (1971) pp. 24-27. 

• A.E. Mayo, Evidence for the adiabatic invariance of the black hole 
horizon area, Phys. Rev. D58 (1998) 104007 [gr-qc/9805047]. 
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11. MECHANICS OF CONTINUOUS SYSTEMS 



Forward: All the formulations of mechanics up to now have dealt with 
systems having a finite number of degrees of freedom or infinitely countable. 
However, many physical systems are continuous; for example, an elastic 
solid in vibrational motion. Every point of such a solid participates in the 
oscillation (s) and the total motion can be described only by specifying the 
coordinates of all points. It is not so difficult to modify the previous formu- 
lations in order to get a formalism that works for continuous media. The 
most direct method consists in approximating the continuum by a collection 
of discrete subunits ('particles') and then study how the equations of motion 
change when one goes to the continuous limit. 

CONTENTS: 

11.1 Lagrangian formulation: from the discrete to the continuous case 

11.2 Lagrangian formulation for the continuous systems 

11.3 Hamiltonian formulation and Poisson brackets 

11.4 Noether's theorem 
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11.1 Lagrangian formulation: from discrete to the 
continuous case 



As one of the most simple case for which one can go easily from a discrete 
system to a continuous counterpart we consider an infinitely long elastic 
rod doing longitudinal vibrations, i.e., oscillations of its points along the 
rod axis. A system made of discrete particles that may be considered as 
a discrete approximation of the rod is an infinite chain of material points 
separated by the same distance a and connected through identical massless 
resorts of elastic constant k (see the figure). 




We suppose that the material points can move only along the chain. Thus, 
the discrete system is an extension of the polyatomic lineal molecule pre- 
sented in chapter 6 of Goldstein's textbook. Thus, the equations of motion 
of the chain can be obtained by applying the common techniques used in 
the study of small oscillations. Denoting by rji the displacement of the ith 
particle with respect to its equilibrium position, the kinetic energy is 



1 



where m is the mass of each particle. The corresponding potential energy 
is the sum of the potential energies of each of the resorts: 



1 



(2) 



From the equations and (|2|) we get the Lagrangian of the system 



L = T-V 



m Vi -k {r]i+i - ijiY 



(3) 



which can also be written in the form 



L 



m -2 

— -ka 
a 



a 



(4) 



where a is the equilibrium distance between the points. The Euler-Lagrange 
equations of motion for the rji coordinates read 



m 



Vi -ka 



rji+i - 111 



+ ka 



Vi-i 



0. 



(5) 
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The particular form of L in equation and the corresponding equations 
of motion have been chosen as being the appropriate ones for passing to the 
continuous hmit by means of a — > 0. it is clear that m/a turns into mass per 
lenght unit ^ of the continuous system, but the limiting value of ka is not 
so obvious. We recall that in the case of an elastic rod for which Hooke's 
law holds, the enlargement of the rod per unit of length is proportional to 
the force (tension) according to 

F = Yi, 

where ^ is the enlargement per length unit and Y is the Young modulus. 
On the other hand, the relative enlargement of the length a of a discrete 
system is given by ^ = — rji) /a. The corresponding force required to 
act on the resort will be 



F = k - rji) = ka 



Vi+i - Vi 



and therefore na should correspond to the Young modulus de Young of the 
continuous rod. When one passes from the discrete to the continuous case 
the integer index i for the particular material point of the system turns into 
the continuous position coordinate x; the discrete variable r]i is replaced by 
ri(x). Moreover, the following quantity 

r]i+i - Vi _ vix + a) -rj{x) 
a a 
entering Li goes obviously to the limit 

dr] 
dx ' 

when a — > 0. Finally, the sum over the number of discrete particles turns 
into an integral over x, the length of the rod, and the Lagrangian (Q) takes 
the form 




2J \ \dx. 
In the limit a ^ 0, the last two terms in the equation of motion (|5|) are 

dr] 

dx 

and taking again the limit a — > the last expression defines the second 
derivative of ij. Thus, the equation of motion for the elastic rod will be 

which is the known one-dimensional wave equation of propagational velocity 

(8) 
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The equation (^) is the well-known formula for the sound velocity (the prop- 
agational velocity of the elastic longitudinal waves). 

This simple example is sufficient to illustrate the main features of the tran- 
sition from a discrete to a continuous system. The most important fact that 
we have to understand is the role of the coordinate x. It is not a gener- 
alized coordinate; it is merely a continuous index substituting the discrete 
one i. To any value of x corresponds a generalized coordinate r]{x). Since 
ry also depends on the continuous variable i, maybe we have to write it in 
a more precise way than r/ {x.t) by showing that x, equally to t, can be 
viewed as a parameter of the Lagrangian. If the system would have been 
3-dimensional and not one-dimensional, the generalized coordinates would 
have to be distinguished by three continuous indices x,y, z and write in the 
form ri{x,y,z,t). We note that the quantities x,y,z,t are totally indepen- 
dent and show up only in r] as explicit variables. The derivatives of rj with 
respect to any of them could thus be written as total derivatives without 
any ambiguity. The equation (^) also shows the Lagrangian is an integral 
over the continuous index x; in the three-dimensional case, the Lagrangian 
would write 

L = J J J Cdxdydz, (9) 

where C is called the Lagrangian density. In the case of the longitudinal 
vibrations of the continuous rod, the Lagrangian density is 




and corresponds to the continuous limit of Li appearing in equation (^). 
Thus, it is more the Lagrangian density than the Lagrangian that we used 
to describe the motion in this case. 



11.2 Lagrangian formulation for continuous systems 

We note that in eq. (^) the C for the elastic rod dependes on ?/= dr]/dt, the 
spatial derivative of ??, drj/dx; x and t play a role similar to its parameters. 
If besides the interactions between first neighbours there would have been 
some local forces, C would have been a function of rj. In general, for any 
continuous system £ can be an explicit function of x and t. Therefore, the 
Lagrangian density for any one-dimensional continuous system should be of 
the form 

C = cU^,'^,x,t]. (11) 



dx dt 

The total Lagrangian, following eq.(|9|), will be 

Cdx, 
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and Hamilton's principle in the limit of the continuous system takes the 
form ^ 

51 = 5 J Cdxdt = 0. (12) 

From Hamilton's principle for the continuous system one would expect to 
get the continuous limit of the equations of motion. For this (as in section 
2-2 of Goldstein) we can use a varied path for a convenient integration, by 
choosing r] from a family of functions of rj depending on a parameter as 
follows 

7/ (x, t;a) = T] (x, t; 0) + aC, (x, t) . (13) 

where 77 (x, t; 0) is the correct function satisfying Hamilton's principle and 
Q is an arbitrary function of 'good' behaviour that is zero in the extreme 
points of t and x. If we consider I as a function of a, then in order to be an 
extremal solution for the derivative of / with respect to a it should become 
zero in a = 0. By directly deriving / we get 

(H _ f^^ f dCdji d^d_(dq\ ( dC\ d ( dr]\ ^ \ 

(14) 

Since the variation of ry, a^, should be zero at the end points, by integrating 
by parts on x and t we obtain the relationships 



*2 dC d fdri 



t, da V dt 



dt 




and 



^2 d£ d fdr]\ , 
— dx 



d^da \dx 




From this, Hamilton's principle could be written as follows 

^^^dxdti^-^(^\-i-{^\}{^] =0. (15) 




Jx, \ dTJ dt 

Since the varied path is arbitrary the expression in the curly brackets is zero: 

(16) 




This equation is precisely the right equation of motion as given by Hamilton's 
principle. 

In the particular case of longitudinal vibrations along an elastic rod, the 
form of the Lagrangian density given by equation ( [ToD shows that 

dC _ dr] dC _ ^drj _ q 

at dx 
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Thus, as we would have hked, the Euler-Lagrange equation ( |lq ) can be 
reduced to the equation of motion (^). 

The Lagrange formulation that we presented up to now is valid for contin- 
uous systems. It can be easily generalized to two, three, and more dimen- 
sions. It is convenient to think of a four-dimensional space of coordinates 
Xo = t,Xi= X, X2 = y,X3 = z. 
In addition, we introduce the following notation 



dr] 

dxi 



(17) 



Employing this notation and the four x coordinates the Lagrangian density 
( [llD takes the form: 

^ = ^{Vp,Vp,iy,^i^) ■ (18) 

Thus, the total Lagrangian is an integral extended to all three-dimensional 
space: 



L 



C {dxi 



(19) 



In the case of Hamilton's principle the integral is extended to a region of 
the four-dimensional space 



51 = 5 C {dx 



0, 



(20) 



where the variation of the r]p nullify on the surface S entailing the integration 
region. The symbolic calculation needed to obtain the corresponding Euler- 
Lagrange equations of motion is similar to the previous symbolic exercise. 
Let a set of variational functions be 

r]p [xu] a) = r]p (x^) + aC {x^) ■ 

They depend on a single parameter and reduce to r]p (xu) when the param- 
eter a goes to zero. The variation of I is equivalent to put to zero the 
derivative of / with respect to a, i.e.: 



dl 
da 



dC drip ^ dC dr]p^y 



drjp da drjp^i, da 

Integrating by parts equation (|2l]), we get 

dl 

da 



{dxp 



0. 
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^-^(^W-^idx )+f(dx 
dr]p dxy \drip^uj\ da ^ J ^ dxy \dr]p^^ 



'p,V 



da 



and taking the limit a — > the previous expression turns into: 



dl_\ 
da) 







{dXf^ 



d 



dC 



dC 

dr]p dxv [drjp^u 



/dijp\ 
\da Jq 



0. 



(22) 
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Since the variations of each rjp is arbitrary and independent equation (22) 
is zero when each term in the brackets is zero separately: 

d f d£ \ dC _^ 
dxu \dr]p^y j dr]p 

The equations (p^) are a system of partial differential equations for the field 
quantities, with as many equations as p's are. 

Example: Given the Lagrangian density of an acoustical field 

£ = ^ (^^0 V +2PoV -ff-jPoiV- ff)' 

fiQ is the equilibrium mass density and Pq is the equilibrium pressure of the 
gas. The first term of C is the kinetic energy density, while the rest of the 
terms represent the change in the potential energy of the gas per volume unit 
due to the work done on the gas o por el curso de las contracciones y expan- 
siones que son la marca de las vibraciones aciisticas, 7 es el cociente entre 
los calores molares a presion y a volumen constante obtener las ecuaciones 
de movimiento. 

Solution: 

In the four-dimensional notation, the form of the Lagrangian density is 

^ = ^ ifJ-oVifiVifi + 2-Pof?i,j - lPo'ni,i'nj,j) ■ (24) 
From the equation (|23|) the following equations of motion are obtained 



fJ-oVjfiO - lPom,ij = 0, j = 1, 2, 3. (25) 



Coming back to the vectorial notation the equations ( pq ) can be written as 
follows 

Ho—dt^ - jPqVV • r/ = 0. (26) 
ri 

Using the fact that the vibrations are of small amplitude the relative varia- 
tion of the gas density is given by the relationship 

a = —V • f] . 

Applying the divergence and using the previous equation we obtain 
which is a three-dimensional wave equation, where 



is the sound velocity in gases. 
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11.3 Hamiltonian formulation and Poisson brackets 



11.3.1 Hamiltonian formulation 

Hamilton's formulation for continuous systems is similar to that for dis- 
crete systems. To show the procedure we go back to the chain of material 
points we considered at the beginning of the chapter, where for each rji one 
introduces a canonical momentum 

dL dLi 

Pi = — — = a — —■ (27) 
The Hamiltonian of the system will be 

H=piVi-L = a— ^ -L, (28) 

or 

H = a(^V^-L^ . (29) 

Recalling that in the limit a — > 0, L — > £ and the sum in the equation (|2 
turns into an integral the Hamiltonian takes the form: 

^ = /^-(^^-^)- (30) 

The individual canonical momenta pi , given by equation ( p7[ ) , go to zero in 
the continuous limit, nevertheless we can define a momentum density it that 
remains finite: 

Lima^o— = = TT^ • (31) 



The equation ( |30| ) has the form of a spatial integral of the Hamiltonian 
density Ti defined as 

n = Trr] -C . (32) 

Even when one can introduce a Hamiltonian formulation in this direct way 
for classical fields, we should keep in mind the the procedure has to give a 
special treatment to the time variable. This is different of the Lagrangian 
formulation where all the independent variables were considered on the same 
foot. This is why the Hamilton method will be treated in a distinct manner. 
The obvious way to do a three-dimensional generalizeation of the field rjp is 
the following. 

We define a canonical momentum 

K) = If- (33) 
'p 

where tj^ (xj, t) , vr^ {xi,t) together, define the phase space of infinite dimen- 
sions describing the classical field and its time evolution. 
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Similarly to a discrete system we can seek a conservation theorem for vr that 
look like the corresponding canonical momentum of the discrete systems. If a 
given field rjp is a cyclic variable (£ does not present an explicit dependence 
on r/p), the Lagrange field equation has the form of a conservation of a 
current: 

d dC _ ^ 

that is 

^-A|^ = 0. (34) 
at ax. Of] ■ 

Thus, if r]p is cyclic there is a conservative integral quantity 



= y dVTT^ {Xi,t) 



The generalization for the density (eq. (^)) in the case of the Hamiltonian 
density is 

^ {vp, Vp,, , TTp, = r]^ -C, (35) 

where it is assumed that the functional dependence of Vp can be eliminated 
by inverting the eqs. (^). From this definition, one gets 



o^r^^-^if^-jirX^'' '''' 



as a consequence of eq. (|33|). In a similar way, we obtain 
dn _ djl^ a^dv^ d£_ d£_ 

Now, using Lagrange equations, eq. (|37|) turns into 



Due to the occurrence of C, we still do not have a useful form. However, by 
a sim 
have 



a similar procedure as used for getting the terms and for we 



dn dr]^ dc dv^ dc dc 



Thus, by substituting (^9[) in (5|) we get 



dn d f dn 



drjp dxi \drip^i 



(40) 
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The equations ( pq ) y (^Oj) can be rewritten using a notation closer to the 
Hamilton ones for a discrete system. This is possible by employing the 
concept of functional derivative 

A = (41) 

6ip dtp dxi dip^i 
Since Ti is not a function of n . the equations ( ^6|) and ( |40| ) can be written 



V = , TT = . (42) 



Now, by employing the same notation, the Lagrange eqs. (23) can be written 
as follows 

d ( dC\ 5C 

0. (43) 



dt \drij ^Vp 

It is fair to say that the almost unique advantage of the functional derivative 
is the similarity of the formulas with the discrete case. Moreover, one can 
see the parallel treatment of space and time variables. 

11.3.2 Poisson brackets 

We can get other properties of Ti by developing the total time derivative of 
eq. (^), remembering that es funcion de Vp^Vpji'^p Y ^/i- Thus, we have 
that 

dn . . dV^ dC . dC dV, dC dC 
=TT r] +7r ^ r] ^ . 

dt " " " dt Btjp " dVp dt dT]p ^ dt dt 

In this expression, the second and the forth terms nullifie each other because 
of the definition (|33[). The derivative simplifies to 

=7r V r? ^ . (44) 

dt " " d7]p dri^ ,^ dt dt ^ ' 

On the other hand, considering as a function of VpiVpj^'^p ^i^d vr^, the 
total time derivative is 

(£{_. dn dn ■ dn_dr]p^ dn 
~dt~'''d^p^d^p'^'^dv^i^^~dt' ^ ' 

where we wrote the expression in such a manner to get an easy comparison 
with the second terms of eq. (0), and where using the eqs. (|36D, ( |37D and 
(|39|) we obtain 

^ = -w ^''^ 

which is an analog to the corresponding one for discrete systems. 
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On the other hand, the equality of the total and partial time derivatives does 
not hold. Using Hamilton's equations of motion (eq. ( p6| ) yand([40|)) and 
interchanging the order of derivation, the eq.(^) can be written as follows 

<m _ dTi^^ ( drL\ dr(_dri dH 
dt diTp dxi \ drjp^i J drjp^i dxi dt 

Now, employing eq. (^6|) and combining the terms we finally have 

(m__d_(. dn\ dn 

'dt~d^iy'd^i)^~dt' ^ ^ 

which is the closest we can approximate to the corresponding equations for 
discrete systems. 

When C does not depend explicitly on time t, it will not be in 7i as well. 
This implies the existence of a conservative current and consequently the 
conservation of an integral quantity, which in this case is 



H = J HdV . (48) 

Thus, if 7i is not an explicit function of time, the conserved quantity is not 

but the integral H. 
The Hamiltonian is just an example of functions that are volume integrals 
of densities. A general formalism can be provided for the time derivatives of 
such integral quantities. Consider a given density lA and let it be a function 
of the coordinates of the phase space {r]p,TTp), of its spatial gradients and 
possibly on x^: 

U =U {r]p,iT^,r]p^i,TT^^,Xp^ . (49) 
The corresponding integral quantity is 

U{t) = J UdV (50) 

where the volume integral is extended to all the space limited by the contour 
surface on which r/p y vTp take zero values. Doing the time derivative of U 
we have in general, 

• ■ ^ ^ ^\dV (bl) 

dr]p " drjp^i '''' diTp d-Kp^i dt J 




Let us consider the term 

■= fdV-^ — 

drjp^i J dr]p^i dx.. 



Integrating by parts and taking into account the nullity of rjp andd its deriva- 
tives on the contour surface, we have 

dV- ??p,,= - I dVr]^^ 



dr/p^i ' J " dxi \drii 
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For the term in vTp j one uses a similar technique. Substituting the obtained 
expressions and grouping appropriately the coefficients of V and vip, respec- 
tively, and using the functional derivative notation equation ( |5l] ) is reduced 
to 

dU f { 5U . 5U . dU\ 

Finally, introducing the canonical equations of motion (^), we have 

— = dV < —- 1—1— } + / dV— . (53) 

dt J \^ drjp diTp drjp diTp j J at 

The first integral in the rhs corresponds clearly to the Poisson brackets. If 
lA and W are two functions of density, these considerations allows us to take 
as the definition of the Poisson bracket for integral quantities as 

We define the partial time derivative of U through the following expression 



Thus, the eq. (53) podra could be written 



which exactly corresponds, in this notation, to the equation for the discrete 
systems. Since by definition the Poisson bracket of H with itself is zero, the 
eq. (E^) turns into 

^JL-^JL (57) 

dt - dt ' 

which is the integral form of eq. (|47|). Thus, the Poisson bracket formalism 
ocurr as a consequence of the Hamiltonian formulation. However, one cannot 
perform a description in terms of Poisson brackets for field theories by a step 
by step correspondence with the discrete case. 

Nevertheless there is one way to work out the classical fields which includes 
almost all the ingredients of Hamilton's formulation and Poisson brackets of 
the discrete case. The basic idea is to replace the continuous spatial variable 
or the continuous index by a countable discrete index. 
The requirement that r/ be zero at the end points is a contour condition that 
might be achieved physically only by fixing the rod between two rigid walls. 
Then, the amplitud of oscillation can be represented by means of a Fourier 
series: 

n=0 
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Instead of the continuous index x we have the discrete r/. We could use this 
representation of x only when rj [x) is a regular function, which happens for 
many field quantities. 

We assume that there is only one real field r/ that can be developed in a 
three-dimensional Fourier series 

fc=0 

Here, k\s a, wave vector that can take only discrete modulus and directions, 
so that in only one lineal dimension there will be an integer (or sometimes, 
half-integer) wavelengths. We say that k has a discrete spectrum. The scalar 
subindex k represents a certain order of the set of integer subindices that are 
used to enumerate the discrete values of k; V is the volume of the system, 
which appear as a normalization factor. 

The orthogonality of the exponentials in the entire volume can be stated 
through the relationship 



1 r i[ k-k ,.. , ^ 

-je\ J dV = 5,^,, . (60) 

As a matter of fact, the allowed values of k are those for which the condition 
]) is satisfied, and the coefficients (t) are given by 

* = ^ / ^"^'^V t) dV . (61) 
Similarly, for the density of the canonical momentum we have 

^{-^.t) = ^^Y.Vk{t)e-^'^-^ (62) 

k 

with pk (t) defined as 

Pf' = ^ / ^"^'^^ dV . (63) 

Both and are integral quantities. Thus, we can look for their Poisson 
brackets. Since the exponentials do not contain the fields we have, according 

to m 



V 

that is, by equation (|60|), 



1 [ .ni„-it-T 

Srj 6-ir Sr] Sir 

^ I'dVe-'^-^ 



[Qk,Pk'] = Km' ■ (6^) 
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From the definition of the Poisson brackets it is obvious that 

[(lk,qk'] = [Pk,Pk'] = ■ (65) 
The time dependence of is sought starting from 

1 f Sn6r] 



5r] 6ir 5rj 5tt 



that is 

^^^t) = -^le~^^-^'^. (66) 



yi/2 J §t: 

On the other hand, we have 



dpk J dn dpk 

and therefore we get 



d^l^Tr- (67) 



dpk 

Comparando las ecuaciones ( |67|) y (|6^) tenemos 



- ' ^ . (68) 



* « = IS ■ c'^) 

In a similar way we can obtain the equation of motion for pk 

dH , , 

Thus, PkY Qk^ obey the Hamilton equations of motion. 



11.4 Noether's theorem 

We have already seen that the properties of the Lagrangian (or of the Hamil- 
tonian) imply the existence of conservative quantities. Thus, if the La- 
grangian does not contain explicitly a particular displacement coordinate, 
the corresponding caonical momentum is conserved. The absence of an ex- 
plicit dependence on a coordinate means that the Lagrangian is not changed 
by a transformation that alter the value of that coordinate; we say that it 
is invariant or symmetric for that transformation. 

The symmetry under a coordinate transformation refers to the effects of an 
infinitesimal transformation as follows 

where the variation 5x^ can be a function of all the other Xy. Noether's 
theorem deals with the effect of the transformation of the field quantities 
itselves. Such a transformation can be written 

V {Xf,) Vp {x'p) = Tip {x^) + 5Vp {x^) . (72) 



129 



Here Stj^ (x^) is a measure of the effect of the variations of x^ and of tj^. It 
can be a function of all the other fields r]^ . The variation of one of the field 
variables in a particular point of space x is a different quantity 6r] : 



% = rj^ (xj + 67]^ (xj . (73) 

The characterization of the transformations by means of infinitesimal vari- 
ations, starting from untransformed quantities means that we consider only 
continuous transformations. Therefore, the symmetry under inversion of the 
three-dimensional space is not a symmetry of the continuous type to which 
Noether's theorem can be applied. As a consequence of the transforma- 
tions both in the coordinates and the fields, the Lagrangian will, in general, 
appear as a different function of the field and spacetime coordinates: 

^ iVp K) ' K) ' ^ J ^ ^' {% (^1) ' (^1) ' 2^1) • (74) 

The version of Noether's theorem that we shall present is not of the most 
general form possible, but makes easier the proof, without loosing too much 
of its generality and the usefulness of the conclusions. We shall suppose the 
following three conditions: 

1. The spacetime is Euclidean, meaning that the relativity is reduced to 
the Minkowski space, which is complex but flat. 

2. The Lagrangian density is of the same functional form for the trans- 
formed quantities as for the original ones, that is 

^' a (^1) ' i,^ (^m) ' ^1) = ^ a (^m) ' i,.^ (^m) ' • (75) 

3. The value of the action integral is invariant under the transformation 

y (dx,) C (j]'^ (x^) , r]^^ (x'^) ,x^) = j^C {r]^ (x J , rj^^^ (xj , x^ 

(76) 



Combining the equations dTq) and (7t) we get the condition 



, {dxj C {vl (x^) , r/^,, (xj , x^) - J^C (r?^ (xj , i]^,^ (x J , xj = . 

(77) 

From the invariance condition, the equation (|7^) becomes 

,dx^£(^i]' ,x^'^ - dx^C{r],x^) (78) 

i/ i/ S 

Here, C {r],x^) is a shorthand notation for the total functional dependence, 
S is the three-dimensional surface of the region Q and Sx is the difference 
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vector between the points of S and the corresponding points of the trans- 
formed surface S . The last integral can be transformedthrough the theorem 
of four-dimensional divergence. This leads to the following invariance con- 
dition 







^dx l[£(??',x,) -£(7?,xj] (79) 



Now, using the equation (|73|), the term in the brackets can be written in the 
first-order of approximation as follows 

^ a K) 'i,. K) - K) ,^?p,. K) = + ^H..- 

Using the Lagrange field equations 

Then, the invariance condition ( [79|) ocurr as 

which already has the form of an equation for the conservation of a current. 
It is useful to develop more the condition giving the form of the infinites- 
imal transformation as a function of the R infinitesimal parameters £rr = 
1, 2, R, such that the variations of y 77^ be lineal in Sr- 

Sx^ = SrXru, Sr]^ = er"^^^. (81) 

By substituting these conditions in eq. (|80| ) we get 

Since the Sr parameters are arbitrary, there are rconservative currents as 
solutions of the differential conservation theorems: 

The equations (^) are the main conclusion of Noether's theorem, telling 
that if the system has symmetry properties fulfilling the conditions (1) and 
(2) for transformations of the type given by the equations (|8ll), then there 
exist r conserved quantities. 



Further reading 

R.D. Kamien, Poisson bracket formulation of nematic polymer dynamics, 
cond-mat/9906339 (1999) 
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